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Abstract 


The  theory,  design  and  construction  of  a  mi croprocessor 
controlled  gamma-ray  CT  bone  scanner  are  described.  This 
scanner,  with  its  ability  to  determine  the  radiological 
density  of  trabecular  bone  at  appendicular  sites  in  the 
skeleton,  provides  an  order  of  magnitude  increase  in 
sensitivity  to  changes  in  bone  mineralization  over  existing 
techniques  such  as  absorpt i ometry  and  radiographic  analysis. 
This  technological  development  could  be  invaluable  to  the 
clinician  for  diagnosis  and  in  the  evaluation  of  treatments 
of  osteopenia.  The  design  and  implementation  of  both  a 
single  and  mu  1 1 i -detector  rotat ion- trans 1  at  ion  scanning 
system  are  detailed  along  with  future  plans  for  a  fan-beam 


rotational  system. 
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I.  Historical  Aspects 


It  has  long  been  the  goal  of  medical  radiology  to  obtain 
clear  cross-sectional  views  of  the  body.  As  early  as  1921 
(1)  conventional  or  focal  plane  tomography  was  being  used. 
Conventional  tomography  involves  the  movement  of  both  the 
X-ray  tube  and  film  (or  even  the  body  itself)  in  a 
trajectory  (Figure  1)  so  that  only  the  desired  layer  remains 
in  focus  while  all  other  layers  are  blurred.  This  type  of 
tomography  which  is  analogous  to  a  nar row-depth-of - f i e Id  in 
photography  has  the  disadvantage  of  the  blurred  presence  of 
unwanted  fields. 

Computed  tomography  (CT),  which  is  the  latest  evolution 
in  a  series  of  techniques  used  in  image  recons truct i on , 
eliminates  these  unwanted  planes  by  limiting  radiation  to 
the  desired  plane.  CT  is  based  on  the  principle  that  if 
radiological  attenuation  is  measured  along  a  sufficient 
number  of  lines  passing  through  a  given  plane,  then  the 
radiological  densities  throughout  that  plane  can  be 
determi ned . 

The  mathematical  basis  for  computed  tomography  was 
first  described  by  Radon  (2)  in  1917.  Given  that 

1.  f(x,y)  denotes  the  attenuation  coefficient  at  a  point 
(x,y)  in  the  plane,  and 

2.  that  L  is  any  line  in  the  plane: 

then,  the  logarithm  of  attenuation  is  equal  to  the  integral 
of  f  along  L;  P ( L )  =  f f(x,y)  ds 

where  s  indicates  a  length  along  L.  The  problem  then  is  to 
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ONLY  POINTS  ON  THE  FOCAL  PLANE  ARE  CONTINUALLY  REGISTERED 
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BLURRED. 


FIGURE  1.  CONVENTIONAL  TOMOGRAPHY 
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invert  the  equation:  that  is,  given  P(L),  to  determine 
f(x,y).  Assuming  P(L)  is  given  for  all  lines  L  and  that 
f ( x , y )  is  continuous  with  compact  support,  Radon  showed  that 
f ( x , y )  could  be  calculated  for  all  points  in  the  plane.  At 
the  time  of  its  development  this  mathematical  basis  for 
image  reconstruct  ion  was  purely  theoretical  in  nature  and  in 
actual  fact  was  derived  by  Radon  while  he  was  working  with 
equations  dealing  with  gravitational  fields. 

It  was  not  until  1956  that  this  theory  was  made  use  of 
by  Bracewell  (3)  in  the  field  of  radioastronomy.  Bracewell 
was  able  to  identify  regions  on  the  sun  which  emitted 
microwave  radiation.  However,  the  antennas  being  used  could 
only  focus  on  thin  strips  criss-crossing  the  solar  surface. 
In  order  to  map  out  the  solar  regions  of  microwave  activity 
it  was  necessary  to  convert  these  strips  into  a 
representat i ve  plane  of  solar  activity  -  the  exact  problem 
shown  by  Radon  to  be  mathematically  solvable.  Bracewell 
developed  both  an  analytical  technique  (equivalent  to 
Radon's)  and  an  iterative  technique  for  image 
reconstruct  ion . 

The  first  medical  application  appeared  7  years  later 
when  Cormack  (4)  published  his  first  experimental  results. 
Cormack' s  work  involved  development  of  a  mathematical 
technique  for  the  accurate  reconstruct  ion  of  an  image  from 
x-ray  projections  and  applications  based  on  measurements  of 
simple  phantoms.  This  technique  for  reconstruct  ion  involved 
the  expansion  of  P  and  f  into  Fourier  series  and  the 
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subsequent  matching  of  Fourier  coefficients. 

The  first  commercial  scanner  was  the  E . M . I .  Head 
Scanner  described  in  1973  by  G.  Hounsfield  (5).  Having 
realized  that  the  r ay- i ntegr a  1 s  essentially  composed  a 
discrete  rather  than  a  continuous  number  of  equations, 
Hounsfield' s  approach  to  reconstruction  involved  an 
iterative  solution  to  inversion  of  the  equation  defined  by 
Pj  =  the  summation  over  i  of  W;-  f;  (Figure  2) 

where  W is  a  weighting  function  of  the  radiological  density, 
f,'  .  This  weighting  function  represents  the  contribution  of 
the  ith  cell  to  the  jth  ray. 

From  these  early  days  CT  has  evolved  from  the  point 
where  only  scans  of  the  head  could  be  performed  to  today 
where  scans  of  the  entire  body  can  be  made.  A  1977  study  of 
the  U.S.  Department  of  Commerce  (6)  predicted  that  by  1980 
there  would  be  some  2500  CT  scanners  in  the  U.S.  and  annual 
related  expenditures  of  $1.7  billion.  Over  a  dozen  different 
companies  were  involved  in  marketing  CT  machines  in  1976: 
machines  able  to  delineate  differences  in  soft  tissue 
density  as  low  as  0.5  %  (7). 

A.  Measurement  of  Bone  at  the  University  of  Alberta 

The  use  of  CT  for  the  measurement  of  bone  density,  as 
developed  in  Zurich  by  Ruegsegger  and  El sasser  (8),  is  a 
relatively  new  phenomenon.  Since  changes  in  bone 
mineralization  occur  much  more  rapidly  in  trabecular  bone 
than  in  cortical  bone  (9),  detection  of  small  changes  in 
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J  th  RAY 


THE  RECONSTRUCTION  MATRIX  IS  AN  N  BY  N  ARRAY  OF  CELLS. 

THE  ACTUAL  RAYS  ARE  STRIPS  OF  FINITE  WIDTH.  THE  CONTRIBUTION 
OF  THE  I  th  CELL  TO  THE  J  th  RAY  IS  DENOTED  BY  Wij. 


FIGURE  2.  RECONSTRUCTION  MATRIX 
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mineral  content  are  only  possible  if  trabecular  bone  density 
can  be  quantified.  It  is  for  this  reason  that  CT  scans  which 
can  be  used  to  produce  images  of  bone  cross  section  are  so 
va 1 uab 1 e . 

Although  the  principles  involved  in  reconstruction  of 
an  image  of  bone  are  the  same  as  in  commercial  whole  body 
scanners,  there  are  a  sufficient  number  of  differences  to 
merit  a  specially  designed  CT  scanner. 

The  Division  of  Biomedical  Engineering  at  the  U.  of  A. 
has  been  involved  in  the  clinical  determination  of  bone 
mineral  content  using  absorpt i omet r i c  techniques  for  the 
past  8  years  (10).  This  clinical  experience  led  to  the 
development  of  a  CT  system  for  the  measurement  of  bone 
mineral  content  at  the  University  of  Alberta.  Based  upon  an 
original  design  by  Elsasser  (11)  an  eight  detector, 
trans 1  at i on-rotat ion  gamma-ray  CT  system  has  been 
constructed.  With  this  scanner,  measurements  can  be  made  of 
the  distal  forearm  using  Iodine  125  as  a  photon  source. 

This  thesis  is  concerned  primarily  with  the  control 
system  for  the  scanner.  A  mi croprocessor  unit  or  MPU 
(Motorola  6800  D2  Kit)  serves  as  a  system  master  clock, 
providing  control  of  scanner  movement  and  for  data 
acquisition  through  the  laboratory  computer  (HP  2100)  and 
CAMAC  systems.  An  MPU  based  controller  was  chosen  as  it  was 
felt  that  the  timing  specifications  on  the  scanner  would  be 
too  tight  to  be  controllable  by  the  HP  running  under  a 
multi-user  operating  system. 
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II.  Bone  Physiology  and  Pathology 
In  the  diagnosis  of  various  diseases  of  bone  it  is  desirable 
to  measure  changes  in  mineral  content  of  trabecular  (spongy) 
bone.  To  understand  this  requirement  it  is  necessary  to 
examine  the  physiological  and  pathological  character i st i cs 
of  bone.  Once  the  need  for  bone  mineral  density  analysis  has 
been  described  an  examination  of  existing  analysis 
techniques  will  be  presented. 

Bone  consists  of  an  organic  matrix  that  is  greatly 
strengthened  by  deposits  of  calcium  salts  (12-14).  Five 
percent  of  this  organic  matrix  is  a  homogeneous  medium 
called  ground  substance  while  the  other  95  %  is  made  up  of 
collagen  fibers.  It  is  these  collagen  fibers  which  give  bone 
its  tensile  strength  while  the  calcium  salts  provide  its 
compress i ona 1  strength.  The  calcium  salts  are  composed 
principally  of  calcium  and  phosphate  and  are  formed  in  a 
crystal  called  hydroxyapatite.  The  crystals  are  between  200 
and  400  angstroms  in  length  and  between  10  and  30  angstroms 
in  diameter.  Under  varying  nutritional  conditions  the  ratio 
by  weight  of  calcium  to  phosphate  in  these  crystals  can  vary 
from  1.3  to  2.0. 

As  well  as  having  the  obvious  structural  function  the 
skeleton  acts  as  a  reservoir  for  many  of  the  body's  ions, 
notably  calcium  and  phosphate.  Bone  contains  99%  of  the 
total  body  calcium  and  90%  of  the  total  body  phosphorus. 

From  0.4  to  1.0  %  of  this  bone  calcium  is  readily 
exchangeable  and  this  exchange  occurs  at  points  in  the  bone 
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undergoing  active  calcium  resorption  (remodelling)  (15). 

Thus  there  is  a  constant  circulation  of  ions  between  the 
skeletal  reservoir  and  body  fluids  and  tissues.  Ions  not 
contained  in  the  bone  have  a  wide  range  of  physiological 
activities.  For  example;  calcium  not  in  the  bone  affects 
enzyme  activity,  cell  membrane  permeability,  cardiac  rhythm 
and  neuromuscu 1 ar  excitability.  The  concentration  of  calcium 
in  the  blood  is  maintained  between  9  and  11  mg  per  100  ml 
and  it  is  the  buffering  action  of  bone  that  provides  this 
control.  Exchange  of  calcium  with  the  bone  is  controlled 
mainly  by  parathyroid  hormone  from  the  parathyroids  and 
calcitonin  from  the  thyroid.  Another  means  of  control  of 
calcium  ion  concentration  is  absorption  from  the  gut. 

Vitamin  ,  for  example,  is  required  for  calcium  absorption 
from  the  gut . 

Bone  is  made  up  of  3  types  of  cells.  These  cells  are 
scattered  throughout  an  interstitial  matrix  of  fibrous 
protein  and  occupy  3%  of  the  total  volume  of  the  bone. 
Osteocytes  are  the  ordinary  bone  cells  and  comprise  the  bulk 
of  the  cells.  They  are  incapable  of  mitosis  and  hence  must 
be  removed  before  new  bone  can  be  formed.  Osteoblasts  are 
the  cells  which  secrete  the  ground  substance  and  collagen 
onto  the  bone  forming  surface.  This  uncalcified  organic 
matrix  is  known  as  osteoid.  The  osteoblasts  produce  an 
enzyme,  alkaline  phosphotase,  which  is  thought  to  act  by 
splitting  organic  phosphate  compounds  thus  upsetting  the 
local  calcium-phosphate  balance  causing  precipitation  of 
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calcium  salts  in  the  ground  substance.  Osteoclasts  or  bone 
phagocytes  are  the  cells  involved  in  the  removal  of  bone. 
This  occurs  by  both  a  phagocytic  action  and  the  release  of 
enzymes  and  acids  which  cause  the  simultaneous  dissolution 
of  mineral  and  matrix.  Thus  bone  is  continually  being 
deposited  by  osteoblasts  and  being  resorbed  at  sites  where 
osteoclasts  are  active.  Osteoblastic  deposition  is 
stimulated  by  physical  stress  while  sex  hormones  and  calcium 
supplements  decrease  both  bone  resorption  and  formation. 

On  a  macroscopic  level  there  are  two  main  classes  of 
bone:  cortical  or  compact  bone,  which  is  the  harder  outer 
bone  and  trabecular  or  spongy  bone  which  is  the  less  dense 
inner  bone.  Average  compact  bone  consists  of  25  %  by  weight 
matrix  and  75  %  calcium  salts  while  trabecular  bone  is  made 
up  of  a  much  higher  percentage  of  matrix.  Due  to  a  greater 
surface  area  it  is  in  this  less  dense  trabecular  bone 
structure  where  physiological  and  pathological  changes  in 
matrix  structure  or  calcium  salt  deposition  occur  more 
rapidly  and  hence  it  is  desirable  to  measure  the 
distribution  of  density  in  trabecular  bone. 

A.  Osteoporosis  and  Osteomalacia 

Two  common  diseases  of  bone  are  osteomalacia  or  "adult 
rickets"  and  osteoporosis.  Osteomalacia  is  a  lack  of  calcium 
salts  and  is  due  to  a  lack  of  calcium  and  phosphates  in  the 
extracellular  fluid  (16).  This  is  often  caused  by  a  lack  of 
calcium  absorption  from  the  gut  possibly  due  to  a  vitamin  D 
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deficiency.  Kidney  disorders  are  one  common  non-nutri tional 
cause  of  osteomalacia. 

Osteomalacia  involves  an  increase  in  the  proportion  of 
uncalcified  bone  (osteoid)  with  a  cor  respond i ng  increase  in 
both  the  number  and  thickness  of  osteoid  seams.  The 
mineralized  bone  mass  is  usually  decreased  due  to  a  decrease 
in  calcium  and  phosphate.  The  bones  are  soft  and  skeletal 
deformities  may  ensue.  Symptoms  of  osteomalacia  include 
weakness,  bone  pain  and  tenderness  as  well  as  radiologic 
evidence  of  decreased  mineralization  and  pseudof r actures  of 
the  long  bones. 

Osteoporosis  is  abnormal  organic  matrix  formation 
(17,18)  and  may  be  caused  by:  (1)  lack  of  use  of  the  bones, 
(2)  malnutrition  to  the  extent  that  sufficient  protein 
matrix  cannot  be  formed,  (3)  lack  of  vitamin  C,  (4) 
postmenopausal  lack  of  estrogen  secretion,  and  (5)  old  age 
(due  to  dysfuntion  of  various  protein  anabolic  functions). 
Osteoporosis  therefore  involves  a  decrease  in  bone  mass  due 
to  the  absence  of  a  normal  quantity  of  bone  caused  by 
deficient  production  of  bone  matrix.  The  mineral  content  of 
the  bone  however,  is  grossly  normal  and  there  is  no  excess 
of  osteoid.  Thus  there  is  an  imbalance  between  osteoblastic 
bone  deposition  and  osteoclastic  bone  resorption.  Clinically 
there  is  a  decrease  in  the  number  of  and  thickness  of  bony 
trabeculae.  The  cortical  bone  later  becomes  thinned; 
however,  it  is  still  relatively  dense  compared  to  the 
decreased  density  of  the  trabecular  bone.  There  is  an 
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increase  in  susceptibility  to  fractures  from  compression 
forces  (eg.  vertebral  bone  compression  fractures). 

The  density  of  bone  has  been  shown  to  decline  in  all 
people  from  the  age  of  20  onwards,  the  process  accelerating 
after  middle  age  (16,18).  To  this  extent  bone  becomes 
osteoporotic  after  50  years  of  age  to  a  greater  extent  and 
somewhat  earlier  in  women  than  in  men. 

Decreased  bone  mineralization,  which  may  accompany 
osteopenia  (general  bone  disorders),  appears  radiological ly 
as  areas  of  decreased  density.  As  such,  it  is  possible  to 
use  radiographic  techniques  to  diagnose  and  to  evaluate  the 
effectiveness  of  treatment  in  diseases  related  to  changes  in 
bone  mineralization. 

"Since  the  disorders  causing  rickets  and  osteomalacia 
are  generally  of  long  standing  and  must  have  profound 
effects  on  matrix  compostion  or  mineral  metabolism  before 
they  become  apparent  clinically,  therapy  may  be  necessary 
for  months  or  years  before  healing  is  complete."  1  The  same 
may  be  said  for  osteoporosis.  Some  of  the  present  treatments 
include  mineral  supplements,  gonadal  hormones  and  physical 
therapy.  Some  of  these  treatments  may  have  other  effects. 

For  example,  pharmacologic  doses  of  vitamin  D  may  cause 
hypercalcemia.  Treatments  then  are  very  dose  dependant 
and/or  expensive  and  hence  any  technique  which  will  provide 
an  early  and  accurate  means  of  diagnosis  and  measurement  of 
changes  in  bone  density  has  potential  for  both  therapeutic 
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and  financial  benefit.  Since  changes  in  bone  density  first 
appear  in  trabecular  bone,  a  technique  is  required  that  is 
able  to  measure  density  of  trabecular  bone. 

B.  Physical  Evaluation  of  Bone 

This  section  will  examine  some  of  the  techniques 
presently  being  used  for  physical  evaluation  of  bone.  These 
include  both  invasive  and  non-invasive  techniques. 

Bone  Biopsy 

Bone  biopsies  involve  the  physical  removal  of  a  sample 
of  bone  usually  from  the  iliac  crest  (hip  bone).  This  sample 
is  then  analyzed,  usually  by  microscopic  examination  of  the 
number  of  trabeculae  and  osteoid  in  a  given  area  within  the 
sample.  Although  this  technique  does  provide  valuable 
information  the  sample  removal  is  an  invasive  technique  and 
is  traumatic  for  the  patient.  Results  obtained  vary  greatly 
from  site  to  site  and  once  bone  has  been  removed  the  same 
site  cannot  be  re-examined.  This  then  precludes  the  use  of 
bone  biopsies  as  an  accurate  means  of  monitoring  both  the 
course  of  any  diseases  such  as  osteoporosis  and  osteomalacia 
as  well  as  the  effectiveness  of  any  treatments.  Thus  the  use 
of  bone  biopsies  is  severely  limited. 

Radiology  ( X-rays ) 

Prior  to  1963,  diagnosis  of  osteoporosis  and 
osteomalacia  predominantly  involved  the  radiological 
evaluation  of  vertebral  skeleton.  Radiographic  techniques  or 
X-rays  show  total  density  of  a  number  of  planes  superimposed 
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together  on  one  view  and  hence  any  one  plane  of  high  density 
will  overshadow  any  of  the  other  planes.  Thus  it  is  largely 
cortical  bone  that  shows  on  X-rays  and  hence  measurements  of 
density  of  trabecular  bone  are  not  feasible.  Accurate 
diagnosis  of  either  osteoporosis  or  osteomalacia  using 
X-rays  requires  the  diseases  to  be  significantly  advanced 
involving  decreases  in  bone  mass  of  greater  than  25%  (16,18) 
and/or  compression  fractures.  Due  to  the  lack  of  precision 
of  this  technique,  longitudinal  studies  (studies  over  time) 
are  difficult  to  perform  and  it  is  virtually  impossible  to 
evaluate  the  progress  of  disease  and  the  effectiveness  of 
treatments.  Radiographic  distinction  between  osteoporosis 
and  osteomalacia  is  also  difficult  and  the  two  diseases  may 
coexi st . 

Absorpt i ometry 

In  1963,  Cameron  and  Sorenson  (19)  reported  an  improved 
technique,  known  as  absorptiometry,  for  the  j_n  vivo 
measurement  of  bone  mineral  content.  Their  method  involved 
the  use  of  a  nearly  monoenerget i c  collimated  source  of 
photons  mechanically  coupled  and  col  linear  with  a  collimated 
scintillation  detector.  The  scintillation  detector,  with 
appropriate  electronics,  is  able  to  measure  both  the  energy 
of,  and  the  number  of  photons  striking  it  in  an  interval  of 
time.  Therefore,  one  can  measure  both  the  number  and  energy 
of  photons  striking  the  detector  with  and  without  an  object 
placed  between  them.  These  measurements  allow  the  amount  of 
energy  absorbed  by  the  object  to  be  calculated.  In  the  case 
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of  bone  the  number  of  photons  removed  from  the  beam  by  the 
object  is  related  to  the  mass  of  bone  present.  The 
measurement  technique  then  involves  movement  of  the  source 
detector  system  along  a  linear  path  across  the  object  with 
transmission  measurements  taken  at  discrete  intervals 
throughout  the  translation  (Figure  3).  This  results  in  an 
absorption  profile  of  the  object,  the  area  of  which  is  a 
function  of  the  number  of  photons  absorbed  by  the  object. 

Although  this  technique  is  similar  to  a  conventional 
X-ray  in  that  radiological  attenuation  is  measured,  the 
improved  equipment  results  in  marked  improvements  in  both 
the  accuracy  and  precision  of  the  measurement. 

Absorptiometry  could  be  considered  a  digital  form  of 
radiology  in  that  quantitative  rather  than  qualitative 
transmission  measurements  are  made.  Measurement  of  bone 
mineral  content  on  a  cadaver  forearm  (20)  using 
absorpt i ometry  was  found  to  have  an  accuracy  to  within  4-7  % 
of  that  obtained  from  measurements  of  bone  weight. 
Reproducibility  of  measurements  on  a  single  individual  were 
found  to  have  a  coefficient  of  variation  of  2.5  %.  These 
improvements  in  accuracy  and  precision  of  absorptiometry 
versus  conventional  X-rays  were  due  to  elimination  or 
reduction  of  such  causes  of  error  as  film  exposure, 
development  and  calibration;  measurement  errors  due  to 
scattered  radiation  and  the  broad  beam  polychromatic  x-rays 
as  well  as  errors  from  soft  tissue  effects.  Variations  in 
quantities  of  soft  tissue  are  reduced  by  submersion  of  the 
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limb  in  water,  thus  providing  a  constant  pathlength  of 
attenuation.  The  radiological  attenuation  of  fat,  collagen, 
water  and  other  organic  components  of  soft  tissue  are  all 
approximately  equal. 

Loss  of  vertebral  skeletal  bone  has  been  shown  to  be 
reflected  in  measurements  of  radial  bone  mass,  in  particular 
in  the  distal  region.  A  study  was  done  (21,22)  of  169  white 
women,  113  of  whom  had  no  fractures  while  56  had 
radiological  evidence  of  spinal  fractures.  Those  women  with 
no  evidence  of  fractures  had  higher  values  of  radial  bone 
mass  with  a  significance  of  98.8  %  for  cortical  bone  and 
better  than  99.9  %  for  trabecular  bone.  Although  total  bone 
mineral  content  was  found  to  be  a  reliable  determinant  for 
diagnosis  of  decreased  mineralization  it  was  also  found  that 
in  cases  of  very  slight  or  exceptionally  heavy  built 
patients,  evaluation  of  mineral  content  divided  by  bone 
width  provided  a  more  useful  parameter.  Absorpt iometry 
studies  involve  measurement  of  a  single  radiation  absorption 
profile.  This  allows  one  to  evaluate  total  bone  mineral 
content  and  thickness  of  cortical  bone  but  quantitative 
evaluation  of  bone  mineral  density  is  not  possible.  Thus, 
absorptiometry  requires  a  measurable  change  in  the  thickness 
of  cortical  bone  or  total  bone  mineral  content  to  evaluate 
disease  processes  and  this  typically  requires  a  period  of 
several  months  to  years. 

I n  V i vo  Neutron  Activation  Ana  lysis  ( I VNAA ) 

Neutron  activation  analysis  involves  irradiation  of  a 
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volume  by  neutron  bombardment  and  the  subsequent  detection 
of  emitted  radiation  to  determine  the  quantity  of  a 
particular  material  present  (23).  Irradiation  of  calcium  in 
bone  results  in  the  production  of  a  number  of  different 
isotopes  of  calcium.  By  measuring  the  radiation  given  off  by 
one  of  these  isotopes  the  quantity  of  that  isotope  can  be 
calculated.  Using  the  ratio  of  the  quantity  of  that  isotope 
to  total  calcium,  total  bone  calcium  can  be  calculated.  This 
technique  although  providing  valuable  information  is  a 
costly  process  (requiring  14  MeV  neutron  sources)  with 
limited  accuracy  and  precision  (of  the  order  of  5  %). 

Compton  Scatter i nq 

Another  technique  which  has  been  used  for  the  j_n  vivo 
measurement  of  bone  mineral  density  is  the  Compton  scatter 
technique  (24).  This  technique  involves  the  irradiation  of  a 
collimator  defined  volume  and  the  measurement  of  Compton 
scatter  from  a  particular  portion  of  this  volume  by  a 
collimated  detector.  From  these  measurements  it  is  possible 
to  determine  the  average  electron  density  of  the  measured 
volume.  The  major  limitations  of  the  technique  relate  to  the 
fact  that  only  a  relatively  small  volume  of  material  can  be 
studied  and  as  such  repositioning  problems  (in  order  to 
measure  the  same  volume  at  a  future  time)  are  a  major 
concern . 

Radi oi sotope  CT 

Although  CT ,  with  its  ability  to  image  cross  sections, 
would  appear  to  provide  a  technique  for  measurement  of  bone 
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density  there  are  a  number  of  problems  with  standard 
commercial  X-ray  scanners  that  preclude  such  measurements. 
These  problems  are  mainly  manifest  in  the  areas  of  beam 
hardening  (due  to  use  of  a  multi-energy  X-ray  source), 
physical  resolution  (due  to  the  degree  of  collimation  of  the 
source  and  detectors)  and  the  sample  spacing  (which  is 
required  to  cover  an  abdominal  area,  leaving  too  little 
matrix  in  an  area  of  bone  for  sufficient  analyses  to  be 
carried  out).  A  technique  able  to  measure  density  of  both 
trabecular  and  cortical  bone  with  a  precision  of  better  than 
1  %  using  a  radioisotope  CT  scanner  was  first  reported  by 
Ruegsegger  and  Elsasser  (8)  in  1975.  The  principles  behind 
this  special  purpose  CT  scanner  are  the  same  as  for 
commercial  whole  body  scanners;  however,  changes  have  been 
made  in  order  to  facilitate  the  measurement  of  bone  rather 
than  soft  tissue  density.  The  technique  involves  performing 
a  scan  of  the  radius  and  ulna  at  a  distal  site  ,  usually 
1/10  the  distance  from  the  ulna  styloid  to  the  proximal 
point  of  the  ulna.  At  this  site  although  some  50  %  of  the 
bone  area  is  trabecular  it  contributes  only  10  %  to  the 
total  bone  mineral  content.  Given  this  relatively  large  area 
of  trabecular  bone  it  is  possible  to  separately  determine 
the  density  of  both  trabecular  and  cortical  bone  at  this 
site.  The  technique  is  non- i nvas i ve ,  not  traumatic  to  the 
patient  and  allows  repeat  measurements  of  the  same  sample  to 
be  performed. 

It  is  also  hoped  that  CT  may  be  used  to  examine  the 
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bone-cement  and  metal -cement  contact  on  artificial  hips 
(25).  Although  this  problem  is  somewhat  distinct  from  the 
measurement  of  bone  density  (and  will  require  an  X-ray  tube 
photon  source)  many  of  the  problems  encountered  are  similar 
(i.e.  measurement  of  radiological  attenuation  in  high 
density  areas  such  as  bone).  Upon  implantation  it  is 
difficult  to  determine  the  degree  of  bonding  of  artificial 
hips  and  subsequent  loosening  may  result. 

Use  of  a  photon  source  with  energies  higher  than  the 
27.5  KeV  of  Iodine  125,  such  as 

Gadolinium  153  (40  KeV,  lOOKeV),  may  allow  the  CT  scanner  to 
be  used  for  examination  of  bone  gaps  and  degradation  in 
patients  who  are  having  Knee  joints  implanted.  A  higher 
energy  source  is  required  because  of  the  quantity  of  bone 
involved  and  hence  the  high  radiological  attenuation.  At 
present  it  is  often  necessary  to  perform  an  exploratory 
operation  to  determine  the  sizes  of  bone  gaps  and  the  amount 
of  degradation  of  the  bones.  It  is  possible  that  the  CT 
Scanner  may  be  used  to  eliminate  this  exploratory  surgery. 

In  order  to  fully  appreciate  the  potential  uses  of  a  CT 
Scanner  it  is  first  necessary  to  have  a  basic  understanding 
of  what  computed  tomography  is  and  how  it  works.  This  area 
will  be  explored  in  the  next  chapter. 


III.  Principles  of  CT  Scanners 
"CT"  or  computed  tomography  (7,26-30)  is  now  the  preferred 
name  for  " CAT "( computer i zed  (assisted)  axial  tomography).  A 
tomogram  is  literally  a  picture  of  a  slice  or  section  at  a 
given  orientation.  Radiation  in  the  form  of  gamma  or  X-rays 
are  passed  through  (transmission  tomography)  or  originate 
(emission  tomography)  from  the  desired  plane  of  view  without 
entering  other  areas.  The  amount  of  radiation  absorbed  along 
a  given  pathlength  is  measured  and  is  an  indication  of  the 
sum  of  linear  attenuation  coefficients  or  radiological 
densities  along  the  pathlength.  A  linear  attenuation 
coefficient  is  a  measure  of  the  absorption  of  radiation  as 
it  traverses  a  distance  in  an  object.  It  should  be  noted 
that  the  attenuation  coefficient  in  any  object  is  dependent 
on  the  material  composition  and  the  incident  photon  energy. 
This  sum  of  attenuation  coefficients  is  measured  along  a 
sufficient  number  of  pathways  so  that  an  image  depicting  the 
radiological  densities  (or  attenuation  coefficients) 
throughout  the  object  can  be  reconstructed .  Since  this  work 
deals  primarily  with  transmission  tomography  any  further 
references  to  CT  will  be  referring  to  transmission  CT . 

To  this  date  there  have  been  5  generations  of  CT 
scanners.  These  generations  are  identified  by  the  type  of 
mechanical  motion  used  and  the  number  of  detectors  employed. 
The  first  two  generations  of  scanners  (Figures  4  &  5)  are 
known  as  single  detector  and  multiple  detector 
translation-rotation  systems.  Scans  are  performed  by 
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translating  the  source-detector  configuration  past  the 
object  at  which  point  they  are  rotated  about  the  object  and 
the  translation  is  repeated  in  the  opposite  direction. 
Measurements  of  ray- i ntegr a  1 s  or  total  radiological 
attenuation  are  continually  measured  during  translation. 

This  movement  results  in  a  "smoothing”  or  blurring  of  the 
object  as  the  actual  measured  data  is  taken  over  a  finite 
area  of  the  object.  Translations  and  rotations  are  repeated 
until  the  sum  of  the  rotations  meets  or  exceeds  the  180 
degrees  required  for  complete  reconstruction.  Multiple 
detector  systems  are  advantageous  in  that  rotations  between 
translations  can  be  made  through  a  larger  angle  and  hence 
the  number  of  translations  and  rotations  can  be  decreased 
thus  reducing  scan  time. 

The  third  and  fourth  generation  scanners  (Figures  6  & 

7)  are  multiple  detector  rotation-only  systems.  In  the  case 
of  the  third  generation  scanner  the  source  and  detectors 
each  rotate  about  the  same  point  through  the  required 
circle.  In  the  fourth  generation  scanner  the  detectors  form 
a  stationary  circle  of  360  degrees  and  only  the  source  is 
rotated.  This  eliminates  the  need  to  rotate  the  detectors 
and  their  associated  high  voltage  and  signal  cables.  The 
major  advantage  of  the  rotation  only  scanners  versus  the 
rotat ion- trans 1  at  ion  scanners  is  the  reduction  in  scanning 
time.  This  reduction  in  scan  time  is  important  as  it  reduces 
image  artifacts  that  can  be  caused  by  patient  movement. 

The  fifth  generation  scanners  (Figure  8)  employ  both 
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FIGURE  6.  CT  GENERATION  3 
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multiple  sources  and  multiple  detectors  and  may  allow 
ultrafast  scanning  to  produce  dynamic  representations  of 
cardiac  motion  without  special  gating  or  synchronization 
techniques.  Such  scanners  are  a  relatively  new  phenomenon 
and  introduce  additional  problems  in  terms  of  expense  and 
cal i br at i on . 

Radiological  attenuation  along  lines  through  the  plane 
is  often  measured  in  a  parallel  beam  format  associated  with 
the  first  two  generations  of  scanners.  This  means  that 
attenuations  are  measured  along  a  series  of  regularly  spaced 
parallel  lines  taken  at  a  number  of  angular  orientations 
(Figure  9).  Each  series  of  measurements  taken  at  a  given 
angular  orientation  is  referred  to  collectively  as  a  view  or 
profile.  In  the  cases  of  the  third,  fourth  and  fifth 
generation  scanners  data  is  collected  in  a  rotational 
format (  Figure  10).  This  data  can  be  reordered  and 
interpolated  in  order  to  create  the  parallel  beam  format 
thus  permitting  the  same  reconstruct  ion  algorithms  to  be 
used  for  all  generations  of  scanners.  However,  special 
algorithms  have  been  developed  for  the  rotational 
format  (31,32). 

A.  Back- Pro ject ions 

The  first  technique  (7)  used  for  reconstruct  ion 
involved  a  method  now  referred  to  as  Back- Pro ject i on .  2  For 


2  Capitalization  is  used  to  distinguish  between 
Back-Projection  as  a  method  of  image  reconstruct  ion ,  and 
back-projection  as  a  process. 
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illustrative  purposes  three  profiles  or  projections  of  an 
object  have  been  taken  (Figure  11)  and  are  then 
back-projected  to  reconstruct  the  object  in  Figure  12.  Each 
profile  is  made  up  of  a  discrete  number  of  measurements 
taken  at  regularily  spaced  intervals.  Reconstruction 
involves  the  back-projection  of  each  profile  onto  a 
reconstruction  matrix.  This  matrix  reconstruction  may  be 
performed  optically  (for  example  on  a  screen),  using  analog 
electronics  or  digitally  in  a  computer.  The  signal  intensity 
according  to  a  given  ray  sum  is  applied  to  all  points  that 
make  up  that  ray  and  this  is  done  for  all  projections  giving 
an  approximation  to  the  original  object.  As  seen  in  Figure 
12  points  outside  of  the  original  object  may  receive  some  of 
the  back-projection  intensity  resulting  in  what  is  known  as 
the  "star"  artifact.  Points  within  the  object  also  receive 
components  from  neighbouring  points  meaning  that  subtle 
differences  in  density  cannot  be  determined.  For  these 
reasons  Back- Pro ject i on  is  no  longer  being  used  as  a 
technique  for  reconstruction.  However ,  it  does  serve  as  a 
basis  for  understanding  the  reconstruction  techniques 
presently  being  used. 

The  technique  now  used  involves  measurement  of 
individual  profiles  and  their  respective  angular 
orientations  and  the  subsequent  storage  of  this  data  in  a 
computer.  The  data  is  then  processed  using  mathematical 
algorithms  and  the  recons tructed  object  is  displayed  on  a 

different  grey  levels  to  represent  the  varying 
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densities.  Although  the  method  depicted  in  Figures  11  and  12 
merely  sums  the  individual  profiles,  mathematical  techniques 
allow  for  cancellation  of  back-projections  and  hence  the 
object  can  be  reconstructed  more  exactly. 

B.  Reconstruction  Algorithms 

The  two  main  classes  of  algorithms  used  for 
r econs t r uc t i on  are  the  iterative  and  analytical  techniques. 

I terat i ve  Technique 

The  iterative  technique  involves  formulating  an 
arbitrary  image,  calculating  the  profiles  for  that  image  and 
comparing  these  prof i les  to  those  obtained  experimental ly. 

The  starting  image  is  usually  a  blank  screen  (all  ray 
profiles  equal  zero)  or  a  circular  grey  object  of  uniform 
density  (all  ray  profiles  equa 1  a  constant ) .  Accord i ng  to 
the  differences  in  the  measured  profiles  and  those  of  the 
arbitrary  image  the  image  is  modified  and  the  profiles  are 
compared  again.  This  process  is  repeated  until  differences 
between  the  theoretical  and  experimental  profiles  are  within 
acceptable  limits.  There  are  three  general  classes  of  the 
iterative  technique  (7).  In  the  simultaneous  correction 
technique  each  cell  that  contributes  to  a  ray  is  altered  and 
all  projections  are  corrected  simultaneously.  This  technique 
tends  to  overcorrect  and  hence  iterations  oscillate  about 
the  correct  value.  In  the  case  of  the  ray-by-ray  correction, 
corrections  are  made  to  all  points  of  one  ray  at  which  point 

taken  into  account  prior  to  making 


these  corrections  are 


■ 


32 


further  corrections  to  a  different  ray.  This  technique  is 
found  to  work  best  when  large  angles  are  taken  between 
consecutive  projections  to  be  corrected.  The  third  technique 
is  a  point-by-point  correction  in  which  each  point  is 
corrected  for  all  rays  that  pass  through  it  and  all  past 
changes  are  embodied  into  future  changes. 

In  each  of  the  three  cases  the  corrective  mechanism  can 
be  either  additive  or  multiplicative.  Additive,  which  is  the 
method  primarily  used,  refers  to  the  fact  that  correction  is 
divided  among  cells  according  to  their  weighting  factor.  In 
the  multiplicative  method,  correction  is  applied  to  cells 
according  to  their  present  density.  In  this  case  a  grey 
starting  level  is  required.  The  iterative  techniques  can 
require  a  great  number  of  iterations  and  computations  and 
hence  computing  time  may  be  a  limiting  factor. 

Analytical  Technique 

The  Analytical  Technique  can  be  subdivided  into  two 
major  categories:  (1)  two-dimensional  Fourier  reconstruct  ion 
and  (2)  filtered  Back- Pro ject i on .  Two-dimensional  Fourier 
recons t r uct i on  is  based  on  the  fact  that  the  Fourier 
coefficients  of  the  image  are  equal  to  the  Fourier 
coefficients  of  the  projections  at  the  same  angle.  This 
means  that  if  the  Fourier  coefficients  of  the  projections 
are  known  then  the  Fourier  coefficients  of  the  image  are 
also  known  and  hence  by  taking  their  inverse  transform  the 
image  of  the  object  can  be  reconstructed .  Until  the 
introduction  of  the  Fast  Fourier  Transform  (FFT)  and 
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improvements  in  high  speed  computers  had  taken  place,  this 
technique  was  not  practical  to  use  due  to  the  number  and 
complexity  of  the  calculations  involved.  Mathematical 
development  of  this  two-dimensional  Fourier  reconstruct  ion 
technique  is  shown  on  the  following  pages. 

Filtered  Back- Pro ject i on  or  the  convolution  technique 
is  based  on  this  two-dimensional  Fourier  reconstruct  ion . 
Papers  by  Shepp  &  Logan  (33)  and  Ramachandran  <5 
Lakshmi nar ayanan  (34)  have  shown  that  this  Fourier 
recons t ruct i on  may  be  viewed  in  the  spatial  domain  as  the 
sum  of  each  ray-integral  times  a  weighting  function  of  the 
distance  from  the  ray  to  the  point  of  reconstruct  ion .  Thus 
reconstruction  involves  prefiltering  each  ray  measurement  by 
a  weighting  function  and  then  back  projecting  these  filtered 
ray  measurements  to  reconstruct  the  object.  This  allows  the 
image  to  be  reconstructed  as  each  ray  measurement  is  taken 
rather  than  having  to  wait  for  all  measurements  to  be 
completed.  This  method  of  reconstruction  is  accurate,  simple 
and  greatly  reduces  computation  time  over  other  available 
methods.  The  mathematical  formulation  of  the  convolution 
technique  from  Fourier  reconstruction  is  also  shown  on  the 


fol lowing  pages . 
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ANALYTICAL  RECONSTRUCTION  TECHNIQUES 

A.  y 


The  1 ine  L( t,0)  is  defined 
\  by  the  equation 

X  cose  +  y  si n0  =  t 


I o  is  the  number  of 
incident  photons 

I  is  the  number  of  photons  after  passing 
through  the  object  along  line  L ( t ,  0 ) 

In  this  diagram  there  are  three  different  coordinate 
svstems.  The  s-t  coordinates  are  obtained  by  rotation 
of  the  x-y  cartesian  coordinate  system  while  the 
oo-e  polar  coordinate  system  uses  w  as  the  radial 
component  and  the  angle  0  as  the  angular  measure. 

A  line  L(t,0)  (which  is  defined  by  the  equation 
x  c o s 0  +  y  sin0  =  t)  is  chosen  so  as  to  pass  through  the 
object.  This  line  depicts  a  ray  integral  measurement 
as  would  be  made  during  a  CT  scan.  In  the  parallel 
beam  format  a  number  of  parallel  ray-integral  measurements 
(profile)  would  be  made  along  a  line  L(t,0)  at  which  point 


. 
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the  angle  0  would  be  altered  and  another  profile 
measured.  It  is  due  to  this  pattern  of  measurement 
that  the  s-t  coordinate  system  has  been  introduced. 

The  use  of  the  polar  coordinate  system  permits 
convolution  algorithms  to  be  used  for  reconstruction. 

This  aspect  will  be  discussed  later  in  the  paper. 

First,  some  definitions: 

(1)  f(x,y)  is  the  linear  attenuation  coefficient  at 
the  point  ( x  , y )  ; 

(2)  Io  is  the  number  of  incident  photons;  and 

(3)  I  is  the  number  of  photons  after  passing  through  the 
object  along  line  L(t,0). 

The  total  attenuation  of  the  ray  along  line  L(t,0)  is 
defined  to  be  P(t,0)  which  is  equal  to 

- 1 n (  I / 1 o )  =  /  f ( x  ,y )  ds 

L 

where  ds  is  an  incremental  length  along  L(t,0). 

The  one-dimensional  Fourier  transform  of  P(t,0)  is 

A  00  °° 

P(o),0)  =  /  exp(-ioot)  P(t,0)  dt  =  /  /  f(x,y)  exp(-iwt)  dsdt. 

oo  "°° 

Now,  consider  the  transform  from  the  s-t  coordinates 
to  the  x-y  coordinates  where 

x  =  t  cos0  +  s  sin0  ,  y  =  t  sin0  -  s  cos0  , 

t  =  x  cos0  +  y  sin0  and  s  =  x  sin0  -  y  cos0  . 

The  Jacobian  which  arises  from  the  transformation  is 


. 


, 
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3 ( S , t) 

3  s 

Jx 

3s 

3y 

9 

s  i  n0 

-COS0 

3  (  x  ,  y ) 

3 1 

3  t 

COS0 

s  i  n0 

3  x 

3y 

Hence  dsdt  which  is  equal  to  the  magnitude  of  the 
Jacobian  times  dxdy 


3  (  s  ,  t ) 
3  (  x  ,  y ) 


(dxdy)  =  | 1  I  dxdy  =  dx  dy . 


Therefore  P  ( w ,  0  )  may  be  written  as 

a  °° 

PU,e)  =  /  /  f  ( x  ,y)  ex  p  (  -  i  CJ  (x  cose  +  y  sine  ))  dx  dy 

^  CO 

which  coincidentally  is  also  the  two-dimensional  Fourier 

A 

transform  of  x  and  y  (define  as  f(w,0)). 

A  00 

i.e.  f(w,e)  =  // f(x,y)  exp  (-iu  (  x  cose  +  y  sine  ))  dx  dy 


where  w  is  the  spatial  frequency  in  the  direction  of  the 
t  axis. 


Hence,  cos0  =  x/co  and  sin0  -  y/o>  . 

Now,  if  we  consider  the  transforma ti on  from  the  x-y 
cartesian  coordinate  system  to  the  w-0  polar  coordinate 
system  the  Jacobian  arising  from  this  transformation  is 
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1 

9  ( x  ,  y ) 

9_x 

9oo 

9  x 
90 

cose 

-oo  s  i  n  0 

9  (  CO  ;  0  ) 

9y 

sine 

00CO  S0 

1  9oo 

90 

Hence,  dx  dy  which  is  equal  to  9  (  x ,  y )_ 

9  ( oo ,  0  ) 

=  |  oo  |  doo  d  0  . 

This  means  that  f(x,y)  may  be  written  as 


doo  d 0 


f ( x ,y)  =  1 


TT 


oo  A 


4  TT 


/  d0  /  P(oo,0)  exp(  +i oo  (x  cos0  +  y  sin0  ))|w|d 
2  0 


00 


Thus,  f ( x  ,y )  may  be  uniquely  determined  by  a  knowledge 
of  P  ( t , 0 )  for  all  lines  L.  In  reality  P(t,0)  is  known 
only  along  a  finite  number  of  lines  which  may  be  denoted 
by  t  =  t k  =  ka,  k  =  0,  _1,  _2,  ...  ; 

j  =  j tt  ,  j  =  0,1,2,... 
n 

where  a  is  the  parallel  spacing  between  projections 
and  n  is  the  number  of  views. 


. 


Rewriting  the  integrals  as  finite  sums  we  get 

n  n-l  00  k=  00 

/  d9  *=  Z  tt  and  /  dw"^  Z  a. 

o  j  =  0  n  k=-°° 


By  bounding  the  object  to  a  circle  of  diameter  a 

l<=  co  k-  1  /a 

Z  can  be  rewritten  as  Z 
l^= — oo  k-  1  /  a 


where  there  are  2/a  rays  contained  in  each  profile. 
Hence,  f(x,y)  may  be  written  as 


n-1  k=l  /  a  a  - 

f(Xjy)  =  _a_  Z  Z  P(w,e)  exp( -iw)(x  cos0  +  y  sine  ))  |U>| 

4TTn  j=0  k=-l/a 


' 


Consider  the  Fourier  inversion  formula 


7T  °o  A 

f(x,y)  =  J_  /  d0  /  P(w,0)  exp(iw(x  cose  +  y  sine  ))  /Wj  dw 


Substituting  in  t  =  x  cosQ  +  y  sine  the  inner  integral  may  be 
written  as 


oo  A 

Q(t,6)  =  /  P(W,0)  |W)  exp(itot)  dW  . 

-00 


Now  if  f(x,y)  is  a  smooth  function,  then  P(t,6)  will  be  a  smooth 

A 

function  and  hence  P(W,6)  will  be  bandlimited. 

A 

i . e.  for  U)  >  ft  P(W,0)  ->  0 . 

If  there  also  exists  an  even  function 

A  00 

'F(W)  =  /  T(t)  exp(-iwt)  dw 

— oo 

which  equals  |W|  for  u)  <  Q 
then  Q(t,0)  can  be  approximated  by 

CO  A  A 

Q(t,e)  -  /  t(w)  P(w,e)  exp(iwt)  dw 

— oo 

which  according  to  the  Convolution  theorem  is  a  convolution 
in  the  time  domain. 

oo 

Q(t,e)  =  2 7T  /  T(x)  P(t-x)  dT 


TT  o° 

Therefore,  f(x,y)  =  _J_  /  0  /p(t,0)  'f  (t-i)  dx 

n  — 1 OO 

2tt  0 

Once  again  writing  the  integrals  as  finite  summations  we  get 
n-1  1/a 

f(x,y)  =  _a  Z  E  P(tk,0)  X  cos  0j  +  y  sin  0^  -  tk) 

2n  j=0  k=-l  /  a 
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Thus  f(x,y)  can  be  written  as  a  summation  of  the  ray  integrals 
times  a  weighting  function  T  which  is  a  function  of  the 
distance  between  the  line  L(t^,0)  and  the  point  of 
reconstruction  (x,y).  One  commonly  used  weighting 
function  (derived  by  Shepp  and  Logan)  is 

y(t  )  =  4  ;  T(t,  )  =  -4/ ( TTa 2 ( 4k2 -1 ) )  ,  k  =  +1 ,  _2,  ... 

0  K 


with  linear  interpolations  in  the  intervals. 


Allowing  T  to  be  linear  in  the  intervals  is  not  incompatible  with 

A 

the  requirement  that  for  small  W,  T(W)  -|to|  and  greatly  reduces 
computational  time  since  the  weighting  function  can  be 
computed  and  stored  and  values  at  intermediate  points  can 

be  obtained  through  interpolation. 

The  Shepp,  Logan  weighting  function  T(t)  has  the  Fourier  transform 


A 

t  (to)  = 


sin 


Wa 

2 


(  sin  (Wa/2)y^a/2)) 2 


When  to  ->0,  T(w)  -H2  toal  ((Wa/2)f;a/2))2  -  lw|(l)2-  |to| 

la  2  I 


and  thus  the  necessary  requirement  is  met. 
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Compar i son  of  Analytical  and  I terati ve  Algor i thms 

Both  the  iterative  and  analytical  techniques  have 
character i st i c  advantages  and  disadvantages.  It  is  generally 
accepted  (7)  that  analytical  methods  are  faster  than  are 
iterative.  Use  of  the  convolution  technique  also  allows  the 
image  to  be  reconstructed  as  measurements  are  made.  In  both 
cases  there  are  problems  associated  with  reconstruction.  Due 
to  the  digitization  involved  in  analytical  methods  the 
projections  must  be  bandiimited.  This  means  that  an 
overshoot  phenomenon  wi 1 1  occur  at  sharp  edges  such  as 
bone-flesh  interfaces.  Although  this  effect  can  be  partially 
corrected  by  pref i 1  ter i ng ,  any  such  gains  will  be  offset  by 
a  loss  in  spatial  resolution.  Iterative  techniques  are  able 
to  avoid  the  problems  associated  with  sharp  interfaces. 

With  complete  data  the  two  types  of  algorithms  are 
comparable  with  respect  to  the  amount  of  noise  in  the  image. 
However,  in  cases  where  some  of  the  data  is  missing,  the 
iterative  technique  will  prove  to  be  superior  since  it 
merely  "smooths"  the  image  while  the  analytical  techniques 
assume  the  missing  data  is  the  same  as  the  available  data. 

C.  Photon  Attenuation 

The  ability  of  CT  to  distinguish  between  narrowly 
differing  densities  is  dependent  on  the  character i st i cs  of 
photon  attenuation.  As  photon  radiation  passes  through  an 
absorbing  medium  it  is  absorbed  exponentially  according  to 
the  following  relationship  (35,36); 
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I  (S)  =  Io  exp  -(ju(S-So)) 

where  Io  is  the  number  of  photons  at  a  distance  S=So  and  p 
is  the  linear  attenuation  coefficient.  This  absorption  is 
due  to  four  processes:  the  photoelectr ic  process,  the 
Compton  process,  coherent  scattering  and  pair  production. 

The  linear  attenuation  coefficient  (ju )  has  units  of  inverse 
distance  and  is  a  measure  of  the  fraction  of  photons  removed 
per  unit  length  of  absorber,  p  may  be  converted  to  either 
the  mass,  electronic  or  atomic  attenuation  coefficients  by 
d i v i d i ng  by  either  mass  density  p  ( gm / cm  ) ,  elect  r on  dens i t  y 
(electrons/cm  )  or  atomic  dens i ty  px ( atoms/cm* ) .  At 
moderate  energies  (less  than  1.02  IVIev )  photon  attenuation 
predominantly  consists  of  two  contributions:  (1)  Compton 
scattering,  which  is  proportional  to  electron  density  and 
(2)  photoelectric  absorption  which  is  a  function  of  the 
atomic  number  of  the  element  (approximately  Zy).  At  the 
lower  edge  of  the  energy  spectrum  it  is  the  photoelectric 
effect  which  is  predominant  and  as  such,  measured 
attenuation  is  representat i ve  of  the  atomic  numbers  of  the 
elements  contained  in  the  object.  If  either  simultaneous  or 
sequential  scans  at  different  energies  (37)  were  used  it 
would  be  possible  to  obtain  measures  of  both  the  Compton  and 
photoelectric  effects  allowing  for  evaluation  of  both 
electron  density  and  atomic  number.  Hence  an  j_n  vivo 
elemental  identification  is  possible  using  CT  (38). 
Statistical  Considerations 


Radioactive  decay  is  a  statistical  process  that  can  be 
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modelled  by  a  Poisson  distribution  (39-43).  The 
character i stic  of  a  Poisson  process  is  that  the  mean  and 
variance  are  equal .  This  means  that  the  coefficient  of 
variation  (which  is  equal  to  the  standard  deviation  divided 
by  the  mean)  is  equal  to  1/\j(mean).  One  implication  of  this 
is  that  the  higher  the  associated  count  rate  the  lower  the 
coefficient  of  variation  wi 1 1  be.  For  example,  a  count  rate 
of  10,000  has  an  associated  coefficient  of  variation  of 

1/  v/ 1 0,000  =  1  %. 

Thus  the  lower  the  radioactivity  of  the  source  the 
longer  the  required  counting  period  in  order  to  maintain  the 
same  coefficient  of  variation.  In  order  that  the  statistics 
be  maintained  and  further  blurring  of  the  image  does  not 
occur  (due  to  measurement  over  a  finite  area)  the  scanner 
would  have  to  translate  more  slowly  and  hence  the  overall 
time  of  the  scan  would  be  increased.  If  the  scan  takes  too 
long,  artifacts  due  to  patient  movement  become  a  problem. 

Dead  T ime 

Radiation  detectors  and  the  associated  electronics  have 
finite  resolving  times.  After  a  system  records  a  pulse  it  is 
unable  to  respond  to  another  pulse  for  a  brief  period,  known 
as  the  dead  time.  As  a  result  the  number  of  pulses  recorded 
by  the  counting  system  is  less  than  the  actual  number  of 
pulses  that  occurred.  This  means  that  the  measured 
statistics  of  the  process  have  been  altered  due  to  the 
measurement  equipment. 

In  order  to  rectify  this  situation  it  is  necessary  to 
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measure  the  dead  time  of  the  system  and  to  use  this 
measurement  to  correct  the  measured  pulse  count.  For  most 
applications  system  dead  time  is  independent  of  the  count 
rate.  Hence,  once  system  dead  time  has  been  established, 
corrections  can  be  made  for  any  count  rate. 

System  dead  time  (p)  is  determined  by  measuring  (in 

order ) : 

1.  the  background  count  rate  (B), 

2.  the  count  rate  resulting  from  one  source  (nA), 

3.  the  count  rate  resulting  from  two  sources  (sources  A  & 

B )  ( nj ) ,  and 

4.  the  count  rate  resulting  from  source  B  (ng) . 

If  this  order  of  measurement  is  followed  any  errors  due  to 

accurate  repositioning  of  sources  will  be  eliminated.  If  the 
true  count  rate  is  N  and  np  <  0.05  then  the  following 
expression  is  valid  (44)  ; 

N=  n( 1+np) 

In  the  cases  of  the  above  measurements  this  results  in  the 
following  three  expressions; 

Na+  B  =  n^( 1+nAp) 

Ng+  B  =  ng ( 1 +nQ  p ) 

+  Nj  +  B  =  n$  (  1  +nr  p  ) 

Adding  the  first  two  equations  together  and  subtracting  the 
third  equation  from  the  result  leads  to  the  following 
expression  (solving  for  p); 

p  =  <<v  ■vir "  B)/(n/'ivrV  • 

Thus  for  a  measured  count  rate  (n)  the  fraction  of  lost 
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counts  is  n ( p )  and  hence  n(n)(p)  counts  must  be  added  to  n 
to  give  the  true  number  of  counts  N.  This  correction  is 
easily  performed  by  a  computer. 

Beam  Hardening 

Absorption  of  radiation  is  dependent  on  the  energy  of 
the  radiation.  As  a  polyenerget ic  beam  passes  through  an 
object  the  lower  energy  beams  are  more  severely  attenuated 
than  are  the  higher  energy  beams.  This  results  in  the  beam 
being  "hardened"  as  it  is  composed  of  a  greater  proportion 
of  high  energy  beams  the  farther  along  it  is  in  the  object 
(45-48).  Thus  in  order  to  determine  absolute  densities  of  an 
object  the  beam  must  be  either  monoenerget i c  or  else 
non-linear  corrections  must  be  applied.  Correction  for  beam 
hardening  may  be  either  material  selective  or  non  material 
selective.  Material  selective  correction  takes  into  account 
the  type  of  absorption  material  and  the  resultant  hardening. 
This  requires  knowledge  of  the  types  of  absorption  materials 
present  and  the  quantities  involved  and  hence  an  initial 
reconstruction  must  first  be  made.  Then  the  raw  data  must  be 
manipulated  according  to  the  first  reconstruction  to  produce 
the  corrected  image.  Non-material  selective  beam  hardening 
merely  recognizes  that  beam  hardening  has  occurred  and  uses 
a  beam  hardening  curve  to  correct  the  raw  data  before 
recons  t  ruct i on . 

If  one  measures  the  amount  of  beam  hardening  as  a 
function  of  distance  then  these  curves  may  be  used  to 
correct  raw  data.  For  example,  if  a  measurement  of  In  (Io/I) 
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=  A  (Figure  13)  is  made,  the  corrected  value  B,  which  is  the 
value  that  would  be  otained  if  there  was  no  beam  hardening, 
is  easily  obtained  from  the  curve.  This  procedure  can  easily 
be  performed  for  all  ray i ntegr a  1  measurements  prior  to 
reconstruction.  It  is  interesting  to  note  that  the 
correction  required  for  measurements  of  relatively  low 
values  of  In  (Io/I)  is  small  and  hence  inaccurate 
corrections  for  measurements  such  as  through  skin  (as 
contrasted  to  bone)  introduces  negligible  error. 

Non-material  selective  correction  techniques  for  beam 
hardening  are  obviously  not  as  accurate  as  are  selective 
techniques  however  for  objects  consisting  predominantly  of 
one  or  two  absorption  materials  it  does  provide  a  quick  and 
simple  correction  technique. 

D.  Detector  Types 

Primary  considerations  in  regards  to  selection  of  a 
detector  include  resolution,  efficiency,  energy  of 
radiation,  counting  rates,  and  in  some  applications,  size. 
Energy  resolution  of  a  detector  is  its  ability  to 
distinguish  between  different  energy  photons.  This  is  often 
defined  as  a  percentage  or  ratio  of  the  full  width  half 
maximum  divided  by  the  primary  energy  level.  Detector 
efficiency  refers  to  the  percentage  of  rays  at  a  specified 
energy  that  will  interact  with  the  detector  rather  than 
being  simply  passed  through  the  detector. 

Three  common  detector  types  are  gas, 
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GIVEN  A  LINEAR  ATTENUATION  COEFFICIENT  TJ 
LN( Io/I )  =  JUxS  OR  I  =  Io  EXP( -jJxS) 


FIGURE  13.  BEAM  HARDENING  CORRECTION 


48 


solid  state,  all  of  which  work  on  the  principle  of  photons 
exciting  the  detector  material  resulting  in  the  production 
of  ions,  optical  radiation  and  electrons  respectively.  These 
phenomena  are  then  converted  to  electronic  signals,  the 
magnitude  of  which  are  proportional  to  the  photon  energy, 
and  the  frequency  of  which  are  proportional  to  the  number  of 
photons  per  second  striking  the  detector. 

Solid  state  detectors  have  better  resolution  than  do 
the  other  two  types.  They  are  also  much  more  expensive  and 
hence  are  only  used  in  applications  where  resolution  is 
critical . 

Gas  detectors  although  usually  having  worse  resolution 
and  efficiency  than  scintillation  detectors  do  have  the 
advantage  of  generally  being  more  economical  and  smaller  in 
size.  It  is  for  this  reason  that  most  fan-beam  machines  use 
xenon  gas  detectors  (7)  instead  of  scintillation  detectors. 
Scintillation  detectors  however  are  economical,  readily 
available  and  are  easy  to  use.  These  factors  coupled  with 
reasonable  resolution,  efficiency,  and  counting  rates  have 
led  to  the  utilization  of  scintillation  detectors  for  many 
spectroscopic  applications. 

Sci nt i 1 1  at i on  Detector 

A  scintillation  detector  is  composed  of  two  basic  parts 
:  (1)  a  scintillation  crystal  and  (2)  a  photomultiplier 

tube . 

Scinti 1 1  at  ion  Crystal 

The  scintillation  crystal  is  a  transducer  that 
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converts  photons  of  energy  to  light  pulses  of 
proportional  intensity  to  the  energy  of  the  photon. 

There  are  both  in-organic  and  organic  materials  in  use 
for  this  purpose.  Sodium  Iodide  crystals  Nal(Tl)  are  the 
most  commonly  employed  crystals  for  gamma -ray 
spectroscopy.  It  is  desirable  to  have  a  high  ratio  of 
counts  under  the  full  energy  peak  relative  to  the  total 
number  of  counts  and  this  ratio  is  approximately  four 
times  greater  for  Sodium  Iodide  than  for  example  Cesium 
Iodide  crystals  CsI(Tl),  with  the  same  gamma-ray  primary 
energy.  This  ratio  is  also  increased  by  the  use  of 
collimators  which  restrict  radiation  to  the  central  part 
of  the  crystal.  Use  of  a  Thalium  impurity  (0.1  %)  acts 
as  a  wavelength  shifter  causing  the  crystal  to  emit  2  or 
more  low  energy  photons  (visible  light)  insteao  of  a 
single  Ultra  Violet  photon.  UV  photons  are  not  suitable 
since  they  are  absorbed  by  most  materials  even  those 
transparent  to  visible  light.  At  low  energies,  typical 
resolution  for  Nal  crystals  ranges  from  6-8  %  for 
crystals  under  5  cm  diameter  to  8-10  %  for  crystals  from 
5-10  cm  diameter.  Thicker  crystals  provide  more 
efficient  1  ight  co 1  lection  and  increase  crystal 
efficiency  so  that  the  same  crystal  may  be  used  with 
higher  energy  sources.  If  the  crystal  is  too  thin  high 
energy  rays  may  not  interact  wi th  crystal . 
Photomultiplier  Tube 

A  photoemi ss i ve  cathode  detects  the  light  pulses 
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produced  by  the  crystal  and  in  response  produces  primary 
electrons.  These  electrons  are  electrostatically  focused 
and  accelerated  by  means  of  a  voltage  differential  to 
the  first  dynode  in  the  tube.  Upon  arrival  at  the  dynode 
they  possess  enough  energy  so  that  for  every  primary 
electron  more  than  one  secondary  electron  is  produced. 
This  procedure  is  continued  throughout  an  arrangement  of 
dynodes  (typically  10)  until  the  final  stage  which  is 
the  anode.  The  voltage  at  the  anode  following  the 
passage  of  a  charged  particle  through  the  phosphor  of 
the  crystal  is  equal  to  the  charge  of  the  electrons  at 
the  anode  divided  by  the  capacitance  of  the  anode  and  is 
proportional  to  the  intensity  of  the  light  impulse  over 
a  wide  range  of  photomultiplier  gain.  Head  on  type  PMT  s 
have  a  semitransparent  photocathode  which  is  in  contact 
with  the  inside  of  the  glass  window.  Head-on  versus 
side-on  photomultiplier  tubes  are  generally  used  for 
scintillation  counting  due  to  improved  uniformity  and 
collection  efficiency. 

Anode  Load 

The  signal  on  the  anode  has  a  rise  time  which 
corresponds  to  the  fluorescence  decay  time  of  the 
crystal.  This  signal  decays  through  the  anode  load  RL 
(Figure  14)  which  is  the  resistor  between  the  anode  and 
the  high  voltage  supply  and  the  parallel  combination  of 
the  anode  to  ground  capacitance  (typically  10  pf )  and 
the  output  line  capacitance  to  ground  (typically  10  pf ) . 
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The  capacitances  are  virtually  fixed  values  and  are  not 
easily  changed  and  it  is  only  the  anode  load  that  is 
easily  altered.  The  signal  on  the  anode  must  decay  away 
fast  enough  so  that  pulse  pi leup  is  not  a  problem  and 
slow  enough  so  that  pulse  height  is  significantly  larger 
than  existing  noise.  Pulse  pi leup  occurs  when  signals  do 
not  have  enough  time  to  decay  away  and  the  next  signal 
is  superimposed  on  the  partially  decayed  signal. 

E.  Performance  Character i st  i  cs  Of  CT  Scanners 

Prior  to  this  point  it  has  been  the  component  parts  of 
a  CT  system  that  have  been  reviewed.  This  section  examines 
how  some  attributes  of  the  overall  system  are  related  to 
these  distinct  segments. 

Resolution  of  a  CT  system  comprises  two  distinct, 
although  interrelated  types  of  resolution:  spatial 
resolution  and  contrast  resolution.  Spatial  resolution  of  a 
CT  system  refers  to  the  minimum  size  object  that  can  be 
imaged  by  the  system  when  that  object  is  part  of  a  periodic 
structure.  Contrast  resolution  is  that  percent  change  in 
contrast  at  an  interface  that  can  still  be  imaged  given  a 
minimum  spatial  resolution. 

Spatial  Resolution 

Assuming  a  constant  photon  flux  there  are  four  major 
factors  which  affect  the  spatial  resolution  of  a  CT  system 
(49).  These  are  the  width  of  the  detector  element  aperture, 
the  distance  between  reading  or  sampling  points,  the  form  of 
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the  convolution  filter  used,  and  the  element  size  of  the 
display  picture  (pixel  size).  In  order  to  examine  the 
spatial  resolution  of  a  system  it  is  useful  to  employ  a 
concept  Known  as  the  Modulation  Transfer  Function  (MTF). 
Values  of  MTF  vary  from  0  to  1  and  are  an  indication  of  how 
well  the  imaging  system  transfers  the  frequency  components 
of  the  object  into  the  final  image.  An  MTF  value  of  1  (or  0) 
means  that  all  (or  none)  of  the  frequencies  present  in  the 
original  object  exist  in  the  final  image.  Good  spatial 
resolution  then,  would  be  indicated  by  large  values  of  the 
MTF  extending  into  the  high  frequency  range. 

Aperture  Size 

The  effect  of  aperture  size  can  be  seen  by 
examining  the  results  obtained  by  viewing  a  series  of 
repeating  objects  through  different  size  apertures  while 
holding  all  other  variables  constant.  From  Figure  15  we 
see  that  when  the  size  of  the  aperture  is  much  smaller 
than  that  of  the  object  spacing,  the  amplitude  of  the 
output  signal  approaches  that  of  the  input  signal.  Thus 
the  aperture  limited  MTF  or  the  ATF  would  approach 
unity.  As  can  be  seen  from  Figure  16,  if  the  aperture 
size  is  increased,  the  amplitude  of  the  output  signal 
becomes  smaller  than  that  of  the  input  signal  and  the 
value  of  the  ATF  is  decreased.  When  the  size  of  the 
apertures  approach  twice  that  of  the  object  periodicity, 
all  information  is  lost  and  the  ATF  approaches  zero.  The 
net  effect  is:  the  smaller  the  aperture,  given  other 
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parameters  constant,  the  greater  the  spatial  resolution 
of  the  system. 

Sampling  Rate 

According  to  the  Nyquist  theorem,  in  order  to 
accurately  reconstruct  an  object  that  is  part  of  a 
periodic  structure,  the  sampling  rate  must  be  at  least 
twice  that  of  the  highest  frequency  component  in  the 
structure.  The  tradeoff  associated  with  sampling  at  high 
frequencies  is  a  resultant  increase  in  the  noise  of  the 
reconstructed  image  (33).  If  the  image  contains 
frequencies  greater  than  the  Nyquist  frequency,  aliasing 
will  occur  resulting  in  streaks  in  the  image  emanating 
from  such  sources  as  bone  edges  which  contain  a  number 
of  high  frequency  components. 

Convolution  Filter 

It  has  been  found  (49)  that  the  higher  the 
frequency  response  of  the  convolution  fi iter  used  in  the 
reconstruction,  the  greater  the  resultant  spatial 
resolution.  The  problem  is  that  the  higher  the  frequency 
response  of  the  filter,  the  greater  the  resultant  noise 
of  the  reconstructed  image  and  hence  a  tradeoff  is 
requ i red . 

Contrast  Resolut ion 

Another  area  of  importance  to  CT  scanners  is  their 
contrast  resolution.  This  is  measured  by  examining  the 
system  limiting  spatial  resolution  with  different  contrast 
med ia.  The  limiting  factor  of  CT  syst ems  in  regards  to 
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contrast  resolution  is  the  presence  of  noise.  The  signal  to 
noise  ratio  can  be  improved  by  increasing  the  incident 
photon  fluence  rate  (reflux).  However,  increasing 
radioactivity  will  result  in  a  corresponding  increase  in 
patient  dosage.  One  measure  of  contrast  resolution  used  for 
commercial  scanners  (49)  is  obtained  by  multiplying  the  hole 
diameter  (in  mm)  by  the  contrast  (in  percent)  and  the  square 
root  of  the  dose  (in  rads).  Inclusion  of  the  dose  in  this 
calculation  results  in  a  measure  of  contrast  resolution 
independent  of  dose  rate. 

F.  Commercial  CT  Scanners 

As  mentioned  previously  it  is  necessary  to  make 
alterations  to  commercial  CT  scanners  before  they  can  be 
effectively  used  for  measurement  of  bone  mineral  density. 
Having  examined  some  of  the  theoretical  concepts  of  CT 
scanning  an  examination  of  the  differences  between  a 
commercial  CT  scanner  and  one  designed  for  measurement  of 
bone  mineral  density  will  now  be  given. 

Commercial  whole  body  CT  scanners  are  both  large  and 
expensive:  typically  $600,000  -  $1,000,000  and  the  source  of 
radiation  is  virtually  always  an  X  ~  r  ay  tube.  These  scanners 
are  oriented  towards  measurement  of  soft  tissues  such  as  the 
abdomen.  In  these  areas  the  application  of  the  scanner  is 
generally  to  locate  tumors  which  will  be  depicted  as  an  area 
of  density  different  from  that  of  the  surrounding  region. 
Absolute  densities  are  generally  not  of  concern  so  much  as 
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the  ability  to  distinguish  between  areas  of  different 
densities.  Commercial  CT  scanners  are  typically  fourth 
generation  scanners  allowing  scans  to  be  performed  in  under 
five  seconds.  Increased  speed  reduces  the  problem  of 
movement  artifacts  and  means  less  inconvenience  to  the 
patient . 

In  transmission  CT  the  source  of  radiation  may  be 
either  an  X-ray  tube  or  a  radioisotope.  An  X-ray  tube  is 
advantageous  over  a  radioisotope  due  to  a  greater  photon 
flux.  They  do  however  have  the  disadvantages  of  a  high 
initial  cost  ($30,000),  large  physical  size  and  the 
requirement  of  a  large  amount  of  calibration  and  adjustment. 
There  is  also  the  possibility  of  a  higher  dose  rate  of 
radiation  with  an  X-ray  tube  source. 

There  are  also  many  problems  that  actually  preclude  the 
use  of  commercial  CT  scanners  for  measurement  of  bone 
density.  The  major  problem  is  manifest  in  the  area  of  system 
spatial  resolution.  Use  of  a  multi-energy  X-ray  source  means 
that  beam  hardening  may  be  a  major  problem  and  the  fact  that 
such  scanners  are  looking  at  relative  rather  than  absolute 
density  means  that  the  data  required  for  beam  hardening 
correction  is  often  not  readily  available.  This  is  further 
complicated  by  the  fact  that  the  detectors  used  in 
commercial  CT  scanners  are  not  operated  in  a  pulse  counting 
mode  and  hence  they  cannot  provide  energy  discrimination  of 
detected  photons.  Sources  and  detectors  are  usually  not 
finely  enough  collimated  to  provide  the  initial  physical 
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resolution  and  the  size  of  the  reconstruction  matrix  usually 
means  that  the  area  associated  with  an  object  the  size  of  a 
bone  is  too  small  for  an  accurate  determination  of  average 
density  of  the  area.  This  is  due  to  the  fact  that  the 
reconstruct  ion  matrix  of  such  a  system  is  typically  required 
to  cover  an  area  the  size  of  an  abdominal  cross-section. 

Commercial  CT  scanners  could  be  used  for  the 
measurement  of  bone  mineral  density.  However,  they  would 
first  require  recol 1 imat ion  of  the  source  and  detectors  and 
software  would  have  to  be  altered  to  correct  for  beam 
hardening  and  to  increase  the  matrix  size  of  the  area  of 
analysis.  These  changes  would  require  a  significant 
disruption  of  an  expensive  piece  of  equipment:  a  disruption 
not  presently  justified  given  the  alternative  use  of  a 
smaller,  radioisotope  CT  scanner. 


IV.  System  Design  and  Construction 
This  section  will  cover  five  areas  of  design  and 
construction  of  the  scanner.  These  include:  (1)  mechanical 
operation,  (2)  radioactive  source  and  detectors,  (3)  scanner 
control,  (4)  adjustments  required  for  a  mu  1 1 i -detector 
system,  and  (5)  data  processing.  The  initial  system  design 
and  construction  were  for  a  generation  one  CT  system. 
Subsequent  to  the  successful  completion  of  the  generation 
one  system  the  necessary  modifications  to  convert  to  a 
generation  two  system  were  undertaken. 

A.  Mechanical  Operations 

The  major  design  requirement  insofar  as  mechanical 
operation  of  this  device  is  concerned,  is  the  ability  to 
properly  position  the  source-detector  configuration  in  both 
the  longitudinal  and  rotational  axes. 

Stand 

The  stand  (Figure  17)  was  designed  to  contain  the 
controls  such  as  the  mi croprocessor  unit,  the  power 
supplies,  relays,  source  and  detectors  on  one  portable  unit. 
The  ability  to  position  the  scanner  at  different  heights  was 
required  so  as  to  facilitate  measurement  of  both  arms  and 
legs.  This  was  accomplished  by  mounting  the  scanner  plates 
on  a  board  with  a  12  inch  diameter  hole  cut  in  the  center 
which  slides  vertically  on  ball  bushings  and  is  driven  by  a 
hand  operated  ball  screw. 
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FIGURE  17.  CT  SYSTEM 
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Source-Detector  A 1 iqnment  and  Motion 

The  source  and  detectors  are  mounted  on  the  top  and 
bottom  edge  of  one  side  of  a  1/2  inch  thick  aluminum  plate 
with  an  elongated  12  inch  circular  hole  cut  in  the  center 
(Figure  18).  The  geometry  of  the  source  with  respect  to  the 
detectors  is  fixed.  This  first  plate  is  mounted  upon  a 
second  plate  (Figure  19)  of  1/2  inch  thick  aluminum  with  a 
12  inch  diameter  circular  hole  cut  into  it.  The  top  plate  is 
attached  to  the  second  plate  by  means  of  Thomson  bearings 
mounted  on  the  first  plate  being  guided  along  case  hardened 
steel  rails  fastened  to  the  second  plate.  This  arrangement 
allows  the  first  plate  (and  hence  the  source-detector 
arrangement)  to  translate  linearly  with  respect  to  the 
second  plate.  The  second  plate  is  mounted  onto  a  third  plate 
(Figure  20)  which  contains  a  circular  bearing  with  an 
external  gear  drive  which  is  used  to  rotate  the 
source-detector  arrangment  about  the  object.  These  plates 
and  the  rotational  bearing  are  mounted  onto  the  stand's 
board . 

Mechani ca 1  Dr i ves 

The  drive  for  the  translation  is  a  stepping  motor 
mounted  on  the  second  plate.  This  motor  is  used  to  drive  a 
belt  connected  precision  ball  screw  arrangement  which  is 
attached  to  the  first  plate.  Each  revolution  of  the  ball 
sc r ew  results  in  a  2.5  mm  mo t i on  of  the  plate.  Since,  in  the 
half-step  mode,  the  motor  makes  one  revolution  every  400 
steps  the  resolution  of  one  step  is  2.5  mm/400  steps  or 
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FIGURE  18.  PLATE  #1 
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FIGURE  20.  PLATE  #3 
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6.25  um/step.  The  ball  screw  is  of  such  a  high  precision  and 
hence  expense  that  a  number  of  microswitch  controlled  relays 
have  been  installed  which  shut  off  power  to  the  stepping 
motors  should  the  scanner  be  operated  outside  its  normal 
scanning  range.  The  plates  are  rotated  by  a  second  stepping 
motor  which  drives  a  worm  gear  arrangement  that  is 
mechanically  linked  to  a  spur  gear  which  in  turn  drives  the 
external  gear  on  the  circular  bearing.  The  worm  gear  was 
necessary  in  order  to  reduce  gear  backlash  caused  by  the 
stepping  nature  of  the  motors.  The  gearing  ratios  mean  that 
approximately  200  steps  of  0.9  degrees  each  are  required  at 
the  motor  shaft  to  induce  a  1  degree  rotation  of  the  plates. 
All  of  the  gearing  ratios  in  the  system  were  known  exactly 
except  the  gearing  ratio  of  the  worm  gear  arrangement.  This 
ratio  was  determined  by  counting  the  number  of  motor  steps 
required  to  rotate  the  plates  360  degrees  as  measured  by  a 
dial  gage  with  a  resolution  of  0.0005  inches. 

B.  Radioactive  Sources  and  Detectors 

The  underlying  principle  behind  any  tomographic  scanner 
is  the  fact  that  if  attenuation  of  a  photon  beam  can  be 
measured  along  a  sufficient  number  of  pathways  through  an 
object  then  the  radiological  density  at  any  point  within  the 
object  can  be  determined.  The  radioisotope  photon  source  is 
1.5  Curies  of  Iodine  125  adsorbed  to  a  2  mm  diameter  zeolite 
bead  and  encapsulated  in  aluminum.  The  detectors  are  1/2 
inch  diameter,  1/2  inch  thick  Nal(Tl)  crystals  coupled  to  13 
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mm  diameter  photomultiplier  tubes  with  borosilicate  windows 
and  bi alkali  dynodes .  The  detectors  were  chosen 
predominantly  on  the  basis  of  availability  and  ease  of 
operation.  A  radioisotope  was  selected  due  to  financial 
constraints  and  ease  of  operation.  The  half-life  of  Iodine 
125  is  60.2  days  and  the  source  is  replaced  every  2-3 
half-lives  at  a  cost  of  approximately  $  1 400/source .  The 
mechanism  of  radioactive  decay  of  Iodine  125  is  included  in 
the  appendix  but  basically  a  gamma-ray  of  approximately  27.5 
Kev  is  produced  during  the  decay  process.  This  energy,  being 
relatively  low,  has  a  comparatively  large  linear  attenuation 
coefficient  meaning  that  subtle  differences  in  density  are 
discernible.  At  this  energy  the  total  attenuation 
coefficient  for  Nal  is  approximately  23  per  cm.  Hence  in  a 
1/2  inch  thick  Nal(Tl)  crystal  only  2.06  x  10  of  the  photons 
would  be  passed  through  the  crystal.  For  Gadolinium  153  with 
an  upper  energy  of  100  Kev  the  attenuation  coeff icient  for 
Nal  is  6  per  cm.  and  the  number  of  photons  passing  through  a 
1/2  inch  crystal  would  be  4.91  x  10V  of  those  incident  upon 

the  crystal . 

S i qna 1  Process i nq 

From  the  anode  of  the  photomultiplier  the  signal  is  fed 
into  a  preamp  (Figure  21).  Since  a  Nal(Tl)  crystal  is  used 
the  rise  time  of  the  signal  corresponds  to  the  fluorescent 
decay  time  of  Sodium  Iodide  (250  ns).  Using  a  trial  and 
error  technique  an  anode  load  of  220  ohms  was  determined  to 
be  suitable  in  terms  of  pulse  height  and  pulse  pileup.  The 
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PREAMP  CIRCUIT  USED  WITH  8  DETECTORS  AND  NARROW  HV  CABLES 


OUTPUT  SIGNAL  WITH 
125I  SOURCE 


3  MV 


0  V 

-15  MV 
-150  MV 


NOTE:  -15  MV  IS  THE  UPDATING  POINT  OF  THE  LECROY  623  DISCRIMINATOR. 

DECREASING  R]  SHORTENS  PULSE  WIDTH  AND  SLIGHTLY  DECREASES  PULSE  HEIGHT. 

I.E.  IF  R,  IS  CHANGED  FROM  3.3  Kfi  TO  220  fi:  V  -300  MV  TO  -150  MV 

t  2.2  US  TO  0.85  US 
VQS+1 8  MV  TO  +3  MV 

DECREASING  Cc  AND  Cf  WILL  DECREASE  x  BUT  ALSO  WILL  CAUSE  OSCILLATION. 

THEREFORE,  Cc  AND  C^  WERE  CHOSEN  AS  SMALL  AS  POSSIBLE  WITHOUT 
OSCILLATION. 


FIGURE  21.  PREAMP  CIRCUIT 
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preamp  then  shapes  the  signal,  increases  its  amplitude  and 
transforms  the  impedance.  The  preamp  used  i s  a  Le  Croy  TRA 
1000  preamp  and  is  able  to  provide  a  signal  of  sufficient 
amplitude  so  that  a  linear  amplifier  is  not  required.  This 
particular  preamp  was  chosen  because  of  its  low  input  noise 
(30  pA/  sqrt  (Hz)  r.m.s.)  and  its  frequency  response  (20  MHz 
at  1  mV/uA  gain).  The  output  of  the  preamp  goes  into  a  Le 
Croy  Model  623  8-channel  di scr imi nator  which  provides  an 
output  pulse  for  every  input  pulse  over  a  preset  energy. 

Thus  lower  energy  pulses  are  discriminated  against  and  only 
pulses  above  a  certain  energy  are  counted.  From  the 
discriminators  the  output  pulses  are  sent  to  a  CAMAC  50  M  Hz 
Scaler  which  counts  the  number  of  pulses  received  in  a  given 
interval.  This  counting  interval  can  be  set  by  the  HP  2100 
computer  which  is  interfaced  to  the  CAMAC. 

C.  Scanner  Control 

The  drive  motors  are  two  Slo-Syn  model  M092  FC09 
stepping  motors  with  nominal  torque  ratings  of  200  oz 
inches.  A  stepping  motor  translates  an  electrical  pulse  into 
a  precise  mechanical  motion  of  the  shaft.  This  motion  is  in 
fixed,  repeatable  increments  permitting  accurate  positioning 
of  the  motor  shaft.  The  motors  can  move  in  either  0.9  degree 
or  1.8  degree  increments  per  input  pulse.  The  position  of 
the  shaft  is  to  within  3  %  accuracy  and  this  error  is 
noncumu 1  at i ve  from  step  to  step.  These  motors  require  a  DC 
power  supply  and  associated  electronic  drives  which  contain 
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the  required  circuitry  to  convert  pulses  into  the  proper 
switching  sequences  for  the  motor.  The  power  supplies  are  2 
Slo-Syn  MPS3000  24  V  DC  while  the  electronic  drives  are  2 
Slo-Syn  STM103  3000  step  per  second  motor  control lers 
(Figure  22).  These  controllers  have  internal  oscillators 
that  can  be  used  to  drive  the  motors  in  either  direction  and 
either  in  the  half-step  (0.9  degrees  per  step)  or  the 
full-step  (1.8  degrees  per  step)  modes.  They  are  also  able 
to  receive  pulses  from  external  logic  devices  such  as 
mi croprocessor s  and  minicomputers.  In  the  scanner  external 
logic  sources  include  a  Motorola  mi croprocessor  unit  (MPU) 
and  a  monostable  multivibrator  which  is  used  to  provide 

single  stepping  capability. 

The  control  system  for  this  scanner  was  designed 
Keeping  in  mind  available  equipment  and  computing  facilities 
(Figure  23).  A  Motorola  MEK  6800  D2  mi croporocessor  unit 
(MPU)  provides  the  central  clock  for  the  system.  It  is  used 
to  control  the  two  stepping  motors  which  provide  the 
rotational  and  translational  motions,  as  well  as  to  signal 
the  computer  when  to  collect  data  from  the  scalers.  All  MPU 
software  programs  have  been  assembled  and  are  contained  in 

the  appendix. 

Timing  Cons i derat i ons 

Since  the  mi croprocessor  unit  is  used  to  control  the 
stepping  motors  and  hence  scanner  position  it  is  also 
necessary  that  the  MPU  clock  be  used  to  determine  when  data 
from  the  detectors  should  be  collected.  Upon  receipt  of  a 
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STEPPING  MOTOR  AND  TRANSLATOR  CIRCUIT  DIAGRAM 


FIGURE  22 
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FIGURE  23.  SYSTEM  BLOCK  DIAGRAM 
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signal  from  the  MPU  the  HP  must  respond  quickly  to  the 
signal  and  collect  data  since  the  scanner  is  continually 
moving.  Although  it  was  hoped  that  this  could  be 
accomplished  with  the  HP  running  under  its  normal  multi -user 
operating  system  it  was  found  necessary  to  operate  the  HP  in 
a  stand  alone  mode  while  operating  the  CT  system.  This  was 
due  to  the  fact  that  under  the  multi-user  system  the  time 
from  a  data-ready  signal  (from  the  MPU)  until  data  was 
actually  collected  from  the  scalers  was  variable  over  too 
great  a  range.  The  stand-alone  mode  allows  data  to  be 
collected  within  4  microseconds  (2  HP  clock  pulses)  of  a 
signal  from  the  MPU. 

MPU  Control 

The  Motorola  MPU  is  the  basic  D2  kit  with  the  only 
additions  being  that  of  1/2  K  of  Ram  ( M 6 8 1 0 )  and  2  K  of 
EPROM  (Intel  P2708).  Program  development  and  EPROM  loading 
were  done  on  an  AMI  microcomputer  development  center. 

The  choice  of  the  8  bit  Motorola  as  the  control 
processor  was  based  upon: 

1 .  avai labi 1 i ty  of  parts , 

2.  price,  and 

3.  facilities  for  EPROM  programming. 

Although  the  trend  is  toward  16  bit  and  32  bit 
microprocessors  it  was  felt  that  the  8  bit  microprocessor 
was  sufficient  for  this  particular  application  and  this 
combined  with  the  aforementioned  constraints  led  to  the 
selection  of  the  M  6800. 
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To  facilitate  insertion  and  removal  of  the  object  to  be 
scanned  the  scanner  must  begin  and  end  operation  in  a  "park" 
position  which  may  not  coincide  with  the  actual  scanning 
path.  The  sequence  of  events  which  occur  during  a  scan  are 
as  fol lows : 

1.  the  MPU  determines  the  park  to  start  distance  and 
direction  and  sends  the  appropriate  pulses  to  the  motor 
control lers ; 

2.  a  signal  is  sent  to  the  HP  signifying  that  a  scan  is  to 
begin ; 

3.  a  signal  is  returned  from  the  HP  to  the  mi crocomputer 
unit  signifying  that  the  HP  is  ready  and  that  a  scan  may 
begin ; 

4.  a  signal  is  sent  from  the  MPU  to  the  HP  signifying  that 
a  scan  has  been  started  and  that  acceleration  of  the 
translational  motor  has  begun; 

5.  a  second  signal  sent  to  the  HP  from  the  MPU  signifies 
that  acceleration  has  ceased  and  that  translation 
velocity  has  been  reached; 

6.  all  translation  pulses  are  sent  to  the  motor  controllers 
from  the  MPU; 

7.  pulses  are  sent  from  the  MPU  throughout  the  scan  to  the 
HP  signifying  when  to  collect  data  from  the  detectors, 

8.  after  translation  is  complete  the  MPU  decelerates  the 
translational  motor,  determines  the  required  rotation 
and  sends  the  appropriate  pulses  to  the  motor; 
the  pattern  (2-8)  is  repeated  until  the  appropriate 
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number  of  rotations  and  translations  have  been 
per  formed ; 

10.  at  which  point  the  appropriate  pulses  are  sent  to  the 
motors  so  that  the  scanner  is  returned  to  the  park 
posi t ion . 

Acceleration  and  deceleration  of  the  stepping  motors  is 
required  since  motor  torque  decreases  with  increasing  speed. 
Thus  in  order  to  ensure  that  all  pulses  sent  to  the  motor 
result  in  a  step  while  at  the  same  time  operating  at  the 
maximum  allowable  speed,  acceleration  and  deceleration  of 
the  motors  is  needed.  This  is  accomplished  by  the 
mi croprocessor  which  uses  a  subroutine  ASD  (see  appendix)  to 
provide  a  series  of  pulse  trains  at  varying  frequencies. 

MPU  Based  PIA 

The  mi croprocessor  unit  (MPU)  is  interfaced  to  the 
motor  controllers  and  the  HP  by  an  on-board  integrated 
circuit  (M  6820  PIA  -  peripheral  interface  adapter).  The 
PIA  contains  an  8-bit  bidirectional  data  bus  for 
communication  with  the  MPU,  two  bidirectional  8-bit 
buses  (output  registers)  for  interface  to  peripherals, 
two  programmable  control  registers,  two  programmable 
data  direction  registers  and  four 

individually-controlled  interrupt  input  lines;  two  of 
which  are  usable  as  peripheral  control  outputs.  The  MEK 
6800  D2  Kit  contains  two  PIAs,  one  of  which  is  used  to 
interface  the  keyboard.  The  second  PIA  appears  to  the 
MPU  as  four  memory  locations  located  at  $8004  -  $8007 
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(where  "$"  denotes  a  hexadecimal  number).  The  output 
register  (OR)  and  data  direction  register  (DDR)  of  port 
A  (one  8-bit  bus  used  for  interfacing)  are  located  at 
$8004  and  the  control  register  (CR)  is  located  at  $8005 
while  the  correspondi ng  registers  for  port  B  are  located 
at  $8006  and  $8007.  The  contents  of  the  control 
registers  determine  the  operation  mode  of. the  interrupt 
lines  and  whether  an  output  register  or  data  direction 
register  is  to  be  addressed.  This  latter  feature  is 
required  since  the  OR  and  DDR  share  the  same  address 
location.  The  data  direction  register  determines  whether 
the  cor respondi ng  interface  lines  to  the  output  register 
are  inputs  (DDR=0)  or  outputs  (DDR=1). 

To  control  the  scanner  6  lines  of  the  PIA  port  A 
are  used.  Two  lines  are  needed  for  each  of  the  two 
stepping  motors  (bidirectional)  and  two  lines  are  used 
in  a  handshake  mode  between  the  HP  and  the  MPU.  The  MPU 
is  initialized  by  a  subroutine  (PIA)  which  first  sets 
control  register  A  (CRA)  permitting  DDRA  to  be 
addressed.  DDRA  is  then  set  to  $1F  (bit  7  is  the  most 
significant  bit)  so  that  lines  0-4  of  ORA  are  outputs 
while  lines  5-7  of  ORA  are  inputs.  Control  register  A  is 
then  reset  so  that  ORA  can  be  addressed.  The  normal 
state  of  the  output  register  is  set  equal  to  $1F  (only 
bits  0-4  may  be  set  by  the  MPU).  Signals  to  the 
translation  motor  are  sent  along  lines  3  and  4  ($0F  and 
$17)  while  signals  to  the  rotational  motor  are  sent 
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along  lines  0  and  1  ( $  1 D  and  $  1 E ) .  A  signal  to  the  HP  is 
sent  on  line  2  ( $  1 B )  and  a  signal  from  the  HP  is 
received  on  line  5.  Bits  0-4  are  active  on  the  positive 
transition  meaning  that  subsequent  to  the  normal  state 
of  the  output  register  the  desired  signal  must  be  stored 
in  ORA  and  then  the  normal  state  must  again  be  stored  in 
ORA  at  which  point  the  message  is  transmitted.  This 
pattern  was  used  since  it  was  found  necessary  to 
initialize  the  output  register  at  the  start  of  any  scan 
in  order  to  be  certain  of  the  contents  of  the  register. 
The  signal  from  the  HP  is  level  sensitive  and  responds 
to  a  positive  signal  in  bit  5. 

Acceleration  -  Deceleration  Ramp 

When  operating  the  stepping  motor  from  the  internal 
oscillator,  the  translator  module  automatically 
accelerates  and  decelerates  the  motor  to  prevent  the 
motor  from  missing  steps  or  overshooting  when  stopping. 
This  ramping  of  the  pulse  rate  is  also  necessary  when 
operating  the  translator  module  from  pulses  supplied  by 
an  external  source  (in  this  case  the  MPU).  When  using 
the  internal  oscillator,  acceleration  and  deceleration 
times  to  a  high  speed  are  i ndependant 1 y  adjustable  from 
50  milliseconds  to  1  second.  The  acceleration  and 
deceleration  ramps  are  used  during  the  transition  from  a 
preset  although  adjustable  base  speed  to  a  preset 
although  adjustable  high  speed. 

In  order  that  the  number  of  steps  and  time  for 
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acceleration  and  deceleration  be  Kept  to  a  minimum  Known 
quantity  a  suitable  MPU  controlled  ramp  had  to  be  found. 
The  base  speed  of  the  control  module  was  set  to  zero  and 
using  the  internal  oscillator  a  suitable  ramp  for  a  high 
speed  of  1500  Hz  was  obtained  through  trial  and  error. 

It  was  found  that  the  same  ramp  was  suitable  for  use 
with  both  stepping  motors.  The  upper  speed  of  1500  Hz 
was  selected  since  it  had  been  decided  that  the  maximum 
speed  of  the  motors  during  a  rotation  or  a  translation 
would  be  at  approximately  that  frequency.  Once  a 
suitable  ramp  had  been  obtained  these  pulses  were  then 
measured  with  an  Ortec  ratemeter  (which  converted  the 
pulse  frequency  to  a  proportional  voltage):  the  output 
from  which  was  displayed  on  a  storage  oscilloscope.  Thus 
acceleration  via  the  internal  oscillator  produced  a  ramp 
on  the  oscilloscope  representative  of  the  accelerating 
ramp  used.  The  next  step  was  to  mimic  this  ramp  using 
the  MPU  as  the  pulse  source.  A  subroutine  (ASD  - 
acceleration,  speed,  deceleration)  was  developed  that 
produces  twelve  ramps  of  varying  lengths  and 
frequencies.  As  a  test  of  the  MPU  generated  ramp  the  MPU 
output  was  also  measured  with  the  Ortec  ratemeter  and 
displayed  on  the  oscilloscope  and  the  resultant  waveform 
was  compared  with  that  obtained  from  the  internal 
oscillator  of  the  control  module.  Once  a  close 
approximation  had  been  found  the  MPU  based  ramp  was 


tested  on  the  scanner. 
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The  frequencies,  corresponding  ramp  lengths,  and 
type  and  direction  of  motion  are  loaded  by  ASD  from  RAM 
and  must  be  set  in  RAM  prior  to  calling  ASD.  These 
parameters  may  be  set  so  that  ASD  can  be  used  for  only 
acceleration  or  deceleration  or  to  preset  the  scanner 
location.  A  count  of  all  pulses  sent  during  subroutine 
ASD  is  made.  A  number  of  smaller  subroutines  are  used  in 
conjunction  with  ASD.  These  include  various  acceleration 
and  deceleration  subroutines  which  set  the  parameters  in 
RAM  so  that  ASD  may  be  used  to  accelerate  to  or  from  a 
certain  speed.  There  are  three  deceleration  subroutines 
(DECEL,  DECEL32 ,  DECEL5 )  which  are  used  to  decelerate 
the  scanner  from  speeds  of  less  than  1500  Hz,  less  than 
400  Hz  and  less  than  100  Hz  respectively.  The  respective 
lengths  of  these  three  ramps  are  200  steps,  32  steps  and 
5  steps.  The  upper  speed  limit  of  acceleration  is 
variable  up  to  1500  Hz  merely  by  setting  the  ramp  length 
parameter  contained  in  two  RAM  locations.  For  a  normal 
scan  (with  a  top  speed  of  1318  Hz)  the  frequencies  and 
ramps  are  loaded  from  two  subroutines  (F  and  INIT) 
contained  in  EPROM  into  the  RAM  locations.  The 
frequencies  vary  from  20.2  Hz  to  1115  Hz  and  the 
corresponding  ramp  lengths  vary  from  1  to  42  pulses.  The 
same  ramp  is  used  for  acceleration  and  deceleration  with 
two  hundred  steps  being  required  in  both  cases.  The 
number  of  pulses  at  a  maximum  frequency  of  1422  Hz  can 
be  varied  from  0  to  over  4  million  and  hence  ASD  is 


. 
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suitable  for  presetting  the  scanner  position. 

Another  minor  subroutine  R8D  initializes  parameters 
in  RAM  and  then  calls  ASD  resulting  in  a  25.7  degree 
rotation  of  the  scanner.  This  is  the  rotation  required 
between  translations  in  a  normal  scan  using  the 
8-detector  holder. 

T  r ans 1  at i on 

Once  the  motor  has  been  accelerated  up  to  the 
desired  speed  a  subroutine  (TRANSL)  is  called.  This 
subroutine  makes  use  of  a  number  of  parameters 
previously  loaded  into  RAM  which  determine  the  delay 
between  motor  pulses  and  the  number  of  motor  pulses  per 
signal  to  the  HP.  A  signal  to  the  HP  signifies  that  data 
is  to  be  collected  from  the  scalers. 

Star  t 

At  the  start  of  every  translation  of  a  scan  an  MPU 
based  subroutine  (START)  is  called.  This  routine  (which 
uses  the  handshake  arrangement  between  the  MPU  and  the 
HP)  sends  a  signal  to  the  HP  (signifying  that  a 
translation  is  desired)  and  in  return  waits  for  a  signal 
from  the  HP  (prior  to  allowing  the  MPU  program  to 
continue) .  An  additional  signal  is  sent  to  the  HP  prior 
to  the  start  of  the  ramp. 

Program  Location 

The  MPU  contains  two  EPROMs  (1024  x  8  bit)  located 
at  $6000- $63F F  and  $C000-$C3FF.  The  EPROM  located  at 
$C000  contains  all  subroutines  such  as  ASD  and  TRANSL 
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(see  flowchart  in  appendix)  while  the  main  body  of  the 
control  programs  are  located  in  the  EPROM  at  $6000. 

Since  the  subroutines  are  rarely  altered  this  means  that 
most  changes  to  the  system  control  programs  can  be  made 
by  altering  just  one  of  the  EPROMs.  The  RAM  contained  in 
the  MPU  is  static  RAM  (128  x  8  bit)  and  is  located  at 
addresses  $0000  through  $03FF.  The  reason  that  many  of 
the  control  parameters  are  loaded  by  the  subroutines 
from  RAM  is  that  this  theoretically  makes  it  possible  to 
alter  the  control  parameters  from  the  Keyboard.  At 
present  however,  most  of  the  RAM  parameters  are  loaded 
into  RAM  from  the  main  control  programs  (which  are 
stored  in  EPROM)  at  the  same  time  as  the  main  control 
program  is  run.  This  is  due  to  the  fact  that  the  control 
parameters  are  at  present  seldomly  changed. 

Reoroducibi 1 i ty  of  Measurements 

Bone  mineral  changes  are  significant  along  the 
longitudinal  aspect  of  the  extremities.  Consequently,  for  a 
valid  quantification  of  changes  in  time  in  a  given 
individual,  it  is  important  to  be  able  to  reproduce  the  site 
of  the  scan  (50).  To  achieve  this  required  positioning 
accuracy  a  MPU  controlled  scout  scan  system  was  developed 
which  locates  the  desired  measurement  site. 

The  scout  scan  involves  measurement  of  26  single 
profiles  at  different  points  axially  about  the  measurement 
site.  Reconstruction  of  these  profiles  (see  Figure  24) 
yields  an  X-ray  like  image  which  depicts  the  relative  sizes 
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FIGURE  24.  SCOUT  SCAN 
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and  orientations  of  the  bone(s)  (femur  vs.  radius,  ulna). 

This  information  when  compared  to  the  same  information 
obtained  from  previous  scans  of  the  same  patient  permits 
accurate  location  of  the  desired  measurement  site. 

The  anatomical  criteria  used  to  define  the  scan  site 
are  the  distances  from  the  distal  tips  of  the  radius  and 
ulna  axially  along  the  bones.  At  present,  sites  two 
millimeters  to  either  side  of  the  measurement  site  are  also 
being  evaluated.  These  three  sets  of  data  are  then 
interpolated  in  order  to  provide  a  reconstructed  area 
identical  to  those  previously  measured. 

Single  Stepping 

In  order  to  provide  the  desired  accuracy  for  initial 
positioning  of  the  scanner  it  was  necessary  to  provide  some 
means  for  single  stepping  both  the  rotational  and 
translational  drive  motors.  This  was  accomplished  by  use  of 
a  74123  retr i gger ab 1 e  monostable  mu  1 1 i vi br ator  (51)  as  shown 
in  Figure  25. 

HP /M i croprocessor  Inter faci nq 

The  mi croprocessor  unit  both  sends  and  receives 
positive  signals  of  +  5  volts  while  positive  signals  sent 
from  the  HP  are  +12  volts.  In  order  to  sink  the  required 
current  (12  mA )  to  provide  a  low  signal  to  the  HP  the  output 
from  the  P I A  is  tied  into  two  inverters  before  being 
connected  to  the  HP  (Figure  26).  Positive  signals  from  the 
HP  are  converted  to  +5  volts  by  use  of  a  resistor  zener 
diode  divider  circuit  (Figure  27). 
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FIGURE  25.  SINGLE  STEPPING  CIRCUITRY 
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FIGURE  27.  VOLTAGE  DIVIDER 
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D.  Mu  1 1 idetector  System 

Upon  completion  and  subsequent  testing  of  the 
generation  one  system  the  modifications  necessary  to  change 
the  system  to  a  generation  two  configuration  were 
undertaken.  This  chapter  will  describe  the  problems 
encountered  in  making  these  modifications. 

The  principle  of  the  mu  1 1 i detector  system  is  to 
collimate  a  single  radioactive  source  so  that  an  array  of 
photons  is  produced  which  subtend  a  solid  angle  emanating 
from  the  source.  The  detectors  are  also  collimated  so  that 
scattering  of  rays  between  detectors  is  reduced.  Linear 
scans  are  made  in  exactly  the  same  fashion  as  for  the 
single-detector  systems  however  the  angle  of  each  rotation 
is  n  times  greater  where  n  is  the  number  of  detectors  (this 
assumes  the  detectors  are  spaced  at  the  desired  angular 
increments).  Thus  the  number  of  linear  translations  and  the 
number  of  rotations  can  be  greatly  reduced  decreasing  the 
time  required  for  a  scan.  This  reduction  in  scan  time  is 
significant  in  that  errors  caused  by  patient  movement  can  be 
reduced.  In  this  particular  scanner  moving  from  a  generation 
one  system  to  an  ei ght -detector  generation  two  system 
reduced  the  normal  scan  time  from  approximately  5  minutes  to 
approximately  1.5  minutes.  Although  the  order  of  data 
collection  in  a  mu  1 1 i -detector  system  is  different  from  that 
of  a  single  detector  system  the  reconstruct  ion  algorithm 
used  is  exactly  the  same. 
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Phys i ca 1  L i mi  tat i ons 

In  order  for  radiological  attenuation  of  the  object  to 
be  determined  the  detectors  must  be  allowed  to  measure  what 
is  Known  as  an  "open"  or  unattenuated  count  rate  during  a 
scan.  This  means  that  in  the  mu  1 1 i -detector  system,  for  each 
translation,  all  source-detector  alignments  must  be  totally 
clear  of  the  object  at  some  point  in  the  translation.  Given 
a  finite  translational  pathlength  this  means  that  either  one 
or  both  of  the  object  size  and  the  angle  of  the  detector  arc 
must  be  restricted.  In  this  arrangement  it  is  always 
possible  to  scan  larger  objects  by  not  using  the  outer 
detectors  in  an  array  and  effectively  reducing  the  angle  of 
the  detector  array.  This  would  mean  that  the  rotational 
angle  between  translations  would  have  to  be  reduced. 

Another  limitation  is  the  size  of  the  detectors.  The 
fact  that  these  detectors  can  only  be  packed  together  within 
limits  means  that  in  order  to  have  a  desired  angle  between 
detectors  the  source-detector  distance  cannot  be  less  than  a 
certain  minimum.  This  distance  is  of  utmost  importance  since 
the  photon  flux  impinging  on  the  detectors  follows  the 
i n ve r s e  sguare  law  meaning  that  the  further  the 
source-detector  distance  the  worse  the  available  statistics 
given  the  same  radioactive  source  and  counting  interval. 
Other  considerations  in  regards  to  the  detector  geometry 
include  scan  time,  count  rates  and  available  channels  for 
signal  processing.  The  final  result,  an  8  detector  unit 
similar  to  the  one  depicted  in  Figure  28,  is  a  compromise 
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FIGURE  28.  DETECTOR  HOLDER 


between  the  aforementioned  factors. 
Intervals  and  Subi nterva 1 s 
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In  the  case  of  the  mu  1 1 i detector  system  the  translation 
of  the  scanner  is  perpendi cu 1 ar  in  direction  to  at  most  one 
of  the  source-detector  orientations.  Thus,  while  during  a 
translation  one  detector  may  move  a  distance  of  A 
perpendicular  to  its  orientation,  the  other  detectors  will 
move  distances  somewhat  different  than  "A"  perpendi cu 1 ar  to 
their  own  orientations.  The  reconstruction  algorithm 
requires  that  all  the  detector  spacings  have  to  be  the  same 
distance  apart  and  therefore  some  means  of  correction  has  to 
be  applied.  This  is  done  by  measuring  subintervals;  each 
interval  or  ray-integral  measurement  consisting  of  a  number 
of  sub i n ter va 1 s .  Thus  while  the  detector  that  is  aligned 
perpendicular  to  the  direction  of  motion  will  have  intervals 
consisting  of  a  specified  number  of  subintervals  the  other 
detectors'  intervals  will  be  made  up  of  more  subintervals. 
The  number  of  subintervals  per  interval  varies  from  detector 
to  detector  according  to  their  respective  angular  deviation 
from  the  line  perpendicular  to  the  direction  of  translation. 
Interpolation  between  subintervals  is  typically  required  in 
order  to  accurately  define  an  interval.  In  this  particular 
system  8  subintervals  are  measured  for  every  interval  of  the 
detector  aligned  perpendi cu 1 ar ly  to  the  direction  of  motion. 
Thus  while  in  a  normal  scan  some  1024  subintervals  are 
measured  per  detector  per  translation  they  are  later 
compressed  on  the  HP  to  128  intervals  prior  to 
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reconstruction . 

Center i nq 

The  reconstruction  matrix  requires  that  each  profile 
consist  of  equal  amounts  of  data  collected  about  the  center 
of  the  recons t r uct i on  matrix.  This  is  due  to  the  fact  that 
the  convolution  algorithm  weights  all  ray  measurements 
according  to  their  distance  from  the  point  of 
recons t ruct i on .  In  the  case  of  the  mu  1 1 i detector  scanner 
each  detector  has  a  different  starting  point  so  it  is 
impossible  for  all  of  them  to  collect  equal  amounts  of  data 
about  the  center  of  the  scanner.  In  order  to  compensate  for 
this,  offsets  must  be  calculated  and  applied  to  the 
collected  data  from  each  detector. 

This  problem  is  most  easily  seen  by  examining  the  case 
of  a  single  detector  scanner.  Assume  that  the  scanner  moves 
a  distance  D1  from  the  start  of  a  scan  to  an  object  0  before 
completing  the  scan  (Figure  29).  Then,  if  the  scanner  is 
rotated  through  180  degrees  and  started  once  again  the 
distance  from  the  start  of  the  scan  to  the  object  0  is  now 
D2.  In  order  that  equal  amounts  of  data  be  collected  about 
both  sides  of  the  object  the  data  wi 1 1  have  to  be  offset  by 
an  amount  equal  to  ( D 1  - D 2 ) / 2 .  Thus  if  equal  amounts  of  data 
have  been  collected  about  both  sides  of  the  object  D1  is 
equal  to  D2  and  no  offsets  are  required.  These  calculations 
are  aided  by  a  special  microcomputer  program  (Cent)  which  is 
a  high  resolution  measurement.  An  object  of  high  density  and 
small  physical  size  (i.e.  a  wire)  is  placed  approximately  in 
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RECONSTRUCTION 


IF  THE  PROFILE  IS  TO  CONTAIN  EQUAL  AMOUNTS  OF  DATA  ABOUT 
BOTH  SIDES  OF  THE  OBJECT  THE  POSITION  OF  THE  DATA  WILL  HAVE 
TO  BE  OFFSET  BY  (D-,  -  D2)/2. 


FIGURE  29.  CENTERING  CONSIDERATIONS 
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the  middle  of  the  scan  path.  The  scanner  is  then  translated 
at  low  speed  and  data  is  collected  every  two  motor  steps. 

The  scanner  is  then  rotated  through  180  degrees  and  the  scan 
is  repeated.  Looking  at  Figure  29  we  see  that  information  so 
obtained  will  allow  one  to  calculate  the  offsets  for  each 
detector.  This  procedure  must  be  repeated  every  time  the 
mechanics  of  the  scanner  (such  as  the  source)  are  altered. 

With  an  aluminum-plexiglass  model,  centering  errors  of 
100  motor  steps  or  0.625  mm.  per  detector  are  noticeable 
(Figure  30).  In  the  case  of  a  circular  object  the  centering 
artifact  appears  on  the  image  as  two  circles  of  slightly 
different  radii  being  joined. 

E.  Data  Processing 

Data  for  one  translation  is  stored  in  the  HP  2100  and 
during  a  rotation  is  transferred  to  magnetic  tape.  This  is 
then  processed  and  analyzed  on  the  University  Computer 
(Amdahl  470/  V7)  and  hard  copies  of  the  image  are  obtained 
on  the  associated  electrostatic  plotter  printer.  The 
reconstruct  ion  program  uses  the  convolution  technique  with 
the  filter  function  proposed  by  Shepp  and  Logan  to  produce  a 
256  x  256  reconstruction  matrix.  Prior  to  reconstruct  ion 
both  dead  time  and  and  a  non-material  selective  beam 
hardening  correction  are  made  on  the  Amdahl  as  is  an 
interpolation  between  ray  measurements  and  angles.  The 
interpolation  takes  the  128  data  points  per  projection  and 
interpolates  between  them  in  order  to  produce  256  data 
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points  per  projection  and  also  interpolates  between 
projections  to  produce  112  instead  of  the  normal  56 
projections.  Although  these  interpolations  do  increase 
computation  time  and  have  a  smoothing  effect  on  the  image 
the  increased  matrix  size  is  needed  in  order  that  the 
analysis  program  has  enough  pixels  with  which  to  work. 
Interpolation  between  the  angles  is  necessary  to  remove 
angular  artifacts  from  the  image.  It  has  been  shown  (7)  that 
if  there  are  N  ray-integral  measurements  made  per 
projection,  then  for  the  data  to  be  complete  there  should  be 
M  projections,  where  M=  (*r/4 )  N .  Thus  if  there  are  128 
ray- integral  measurements  made  there  should  be  128  (///4)  or 
approximately  100  projections  made  in  order  for  the  data  to 
be  complete.  Hence  if  the  number  of  angles  were  not 
interpolated  to  produce  112  projections  instead  of  56 
projections,  angular  artifacts  would  be  present  in  the 
reconstructed  image.  The  interpolation  between  measurements 
is  to  produce  256  ray- integral  measurements  so  that 
reconstruction  on  a  256  x  256  matrix  can  be  performed 
easily.  It  should  be  noted  that  it  is  not  necessary  to 
reconstruct  on  a  matrix  the  same  size  as  there  are  number  of 
ray- i ntegr a  1 s  per  profile;  however,  reconstruct  ion  is 

facilitated  if  this  is  the  case. 

In  order  to  use  a  Grinnell  CRT  display,  which  is  tied 
into  the  HP  system  it  is  necessary  to  either  reconstruct  the 
matrix  on  the  HP  2100,  which  takes  about  30  minutes  (for  a 
128  x  128  matrix) ,  or  else  to  load  the  Amdahl  reconst ructed 
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matrix  onto  magnetic  tape  and  to  subsequently  load  the 
matrix  into  the  HP . 

Open  Count 

Although  the  open  count  rate  is  measured  for  all 
detectors  during  each  translation  the  statistics  associated 
with  such  counts  are  of  limited  accuracy  due  to  a  short 
counting  period.  For  this  reason  an  open  count  rate  is  made 
at  the  start  of  the  day.  This  count  is  made  over  a  period  of 
several  minutes  (using  the  Scout  Scan  program)  so  that  the 
coefficient  of  variation  associated  with  the  measurement  is 
low.  This  open  count  (Io)  is  then  used  throughout  the  day  to 
calculate  the  ray-integral  measurements  ( 1 n ( I o/ I ) ) . 

Pixel  S i ze  and  Grey  Leve 1 s 

The  pixel  size  of  the  system  is  determined  by  dividing 
the  scan  length  by  the  number  of  pixels  in  one  side  of  the 

reconstruction  matrix.  The  maximum  size  object  that  can  be 

scanned  is  76.8  mm  and  the  size  of  the  reconstruction  matrix 
is  256  x  256.  Hence  the  pixel  size  of  the  system  is 
76.8  mm/256  or  0.3  mm.  Each  pixel  level  in  the  reconstructed 

matrix  is  defined  by  a  two's  complement  16  bit  number  and 

hence  there  are  32,768  possible  positive  pixel  levels. 
Although  the  system  pixel  size  is  0.3  mm  cor respondi ng  to  a 
256  x  256  matrix  the  hard  copy  imaging  system  is  not  able  to 
resolve  points  of  that  size.  As  a  result  a  128  x  128  matrix 
with  32  grey  levels  is  used  with  a  correspondi ng  pixel  size 
of  0.6  mm.  Since  there  are  four  times  as  many  pixels  in  the 
actual  reconstruction  matrix  as  in  the  image  matrix,  four  of 
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the  actual  pixels  are  averaged  together  to  produce  one  of 
the  image  pixels.  The  Grinnell  display  system  is  able  to 
reproduce  matrices  of  size  512  x  512  with  256  grey  levels 
and  hence  the  system  limiting  pixel  size  for  this  display 
mode  (CRT)  is  0.3  mm. 

Ana lysi s  Program 

While  the  image  depicting  the  bones  may  be 
aesthetically  pleasing  it  is  really  the  parameters  contained 
in  subsequent  analysis  programs  that  are  of  clinical  value. 
The  image's  main  function  is  to  ensure  that  the  scan  was 
done  properly  and  that  no  artifacts  due  to  such  factors  as 
patient  movement  have  been  introduced.  Once  an  image  of  the 
bone  has  been  reconstructed  an  analysis  program  is  used  to 
identify  such  parameters  as  area,  total  mineral  content  and 
average  density.  These  paramenters  are  defined  for  total 
bone,  trabecular  bone  and  the  cortical  or  compact  bone 
(Table  1 ) . 

The  analysis  program  first  defines  where  the  cortical 
bone  begins  and  strips  away  the  surrounding  fleshy  material 
(Figure  31).  This  is  done  by  defining  a  cortical  bone 
density  level  (8000  is  normally  used)  and  defining  all 
levels  of  material  outside  of  those  pixels  with  levels 
greater  than  or  equal  to  8000  as  being  zero.  The  radius  and 
ulna  are  then  separated  and  each  bone  is  analyzed 
separately.  Each  bone  is  then  separated  into  a  number  of 
shells  by  essentially  moving  the  cortical  bone  contour  in, 
one  pixel  at  a  time.  Trabecular  bone  was  initially  defined 
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as  that  area  of  bone  which  comprised  the  inner  50  percent  of 
the  total  bone  area.  It  was  felt  that  this  arbitrary 
definition  of  trabecular  bone  may  not  be  applicable  to  all 
patients.  The  criteria  now  used  involves  first  calculating 
the  average  densities  for  all  shells.  Trabecular  bone  is 
then  defined  as  the  inner  area  of  bone  bordered  by  the 
maximum  change  in  adjacent  shell  densities.  The  CT 
attenuation  coefficients  are  converted  to  bone  density 
(g/ccm)  by  dividing  >u  by  the  average  mass  attenuation 
coefficient  for  bone.  Due  to  the  fact  that  a  non_material 
selective  beam  hardening  correction  is  applied  (as 
contrasted  to  a  material  selective  correction  (48))  these 
measures  of  density  are  only  valid  for  longitudinal  studies. 
Cross  -  sect i ona 1  (i.e.  population)  studies  require  absolute 

values  for  bone  density  and  this  necessitates  material 


selective  corrections. 
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V.  System  Evaluation 


The  scanner  was  evaluated  in  terms  of:  (1)  spatial 
resolution,  (2)  precision,  (3)  accuracy  and  (4)  radiological 
dosage.  The  results  of  these  evaluations  are  discussed  in 
this  section  and  a  technique  for  measuring  contrast 
resolution  is  suggested. 

A.  Spatial  Resolution 

Physical  resolution  of  the  system  was  evaluated  using  a 
2  1/2  inch  diameter  circular  plexiglass  model  with  periodic 
structures  of  various  size  holes  drilled  into  it.  These 
holes  are  drilled  so  that  the  space  between  holes  is  equal 
to  the  hole  diameter.  The  model  was  subsequently  scanned  at 
three  different  speeds  and  the  holes  were  left  open  to  the 
air.  Motor  pulses  during  the  three  speeds  of  translations 
were  sent  at  1318  Hz,  330  Hz  and  82.4  Hz  and  data  from  the 
CAMAC  scalers  was  collected  every  12  motor  pulses  in  all 
three  instances.  Since  these  speeds  are  ratios  of  16:4:1 
respectively  the  resultant  subinterval  counts  obtained 
during  the  scan  were  in  the  ratios  of  1:4:16  respectively. 
Since  coefficient  of  variation  (C.V.)  is  inversely 
proportional  to  the  square  root  of  the  number  of  counts  this 
means  that  the  C.V.s  associated  with  each  of  the 
ray- integral  measurements  were  improved  by  a  factor  of  two 
with  each  decreasing  speed  of  scan.  The  open  count  rates  of 
the  source  for  detectors  1  through  8  were;  on  a  subinterval 
basis:  (1)  374.42,  (2)  349.92,  (3)  334.08,  (4)  345.49,  (5) 
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389.56,  (6)369.65,  (7)  405.93,  and  (8)  379.45.  Thus  the  C.V. 
associated  with  an  average  8  subinterval  per  interval  open 
count  measurement  is  1/  I ( 2948 . 49 )  =  1.84  %.  At  the  medium 
speed  (interval  C.V.  =  0.92  %)  the  spatial  resolution  has 
been  improved  whereas  a  further  decrease  in  the  interval 
C.V.  to  0.46  %  (low  speed)  had  no  apparent  impact  on  the 
spatial  resolution  of  the  system. 

It  should  be  noted  that  one  limit  of  the  system  in 
regards  to  spatial  resolution  is  the  spatial  resolution  of 
the  imaging  system.  In  this  case  although  the  reconstructed 
matrix  is  of  size  256  x  256  with  a  0.3  mm  pixel  size,  the 
hard  copy  image  matrix  is  only  of  size  128  x  128  with  a 
cor r espond i ng  pixel  size  of  0.6  mm.  To  further  complicate 
matters  the  image  display  has  a  magnification  effect  of  1.6 
times  that  of  the  object.  Using  the  Grinnell  display  system 
one  can  display  a  matrix  size  of  512  x  512  and  hence  it  is 
possible  to  view  the  system  limited  resolution  on  the  video 
screen  and  yet  essentially  the  same  results  were  obtained  as 
with  the  hard  copy  display.  At  high  speed  (Figure  32)  the 
system  spatial  resolution  was  limited  to  1/8  of  an  inch 
while  at  medium  speed  (Figure  33)  the  spatial  resolution  is 
improved  to  one  sixteenth  of  an  inch  (approximately  1.5  mm). 
However,  a  further  decrease  in  speed  to  low  speed  (Figure 
34)  had  little  or  no  effect  on  the  spatial  resolution.  This 
is  evidence  that  beyond  a  certain  level  statistics  are  no 

factor  in  regards  to  spatial  resolution. 


longer  the  limiting 


MEASUREMENT 


319 


17- JAN-80 


!■•■■■  ■  ■  ■  ■  I 


■  ■■■■! 


;:s 


-  •  a  I  I  •  •  I  ■■  ■■  I  ■■  II  ■  ■  I  I  I  ■  I  I  I  •  l  •  •  • 

• • ••••■■ . 

■  •  •  ■  ■  BUM*  BH  ■■■!!■  ■■BI3SB  ■■■■■■■■■■■  ■  •  a  •  • 

•  •  •  ■iBaaBaaBDttDKisaaBaBBB9GBiiBBBaB3aaigBiB«i ... 

■  ■  ••■•iiaBlBaaiaiiisBaiBiiiiBaiiiiiaasaaaaiBaaaiiM  • 

• • ■■ aaa ai lalaai Biii<Kai9lliliaaaaaiB9ilaa aaaa a ■ aaa •  « • 

.  . . . . 

•  •  •  ■BaaaiaBBBMaaviiBaBBBVEtiEiuiflaBRCitiaBBaBBaBi.ilaBaBBiaBaaBB  •»  •  •  ■ 

•  ••  . . .  • 

•  •  isaaflBBS3''n"Ba bgobbiibbiii iBSam i ■■■■■■  ■  a  ■  ■  ■  ■  a  a  a  a  8  b  .j  a  a  as ubb a  a  •  •  - 
■  •  •  bBBBB  a  a  axis Jia  aiaaiBBUlai  a  •BfllniMitniMMBi  ■amiBicsaEKSUBai.  • 

> ■laaaBiBBBiasiaiiaiBlfiBlIlBii a  a a ■ ■  ■ ■ a  aaBBB BIS  Ida BBB|aaaaaB aaaa Baaa a 

.  . . .  aaaaBBBQBBBBaaaCSHBSDO  la  a  a  a  a  a  a  a  aaar 

.  . . ■•■•.lam - 

•  •  . . " 

•  . . .  •  ■  a  naCUBaBn - -  - - 

■  BflaBBaiMiiiiMDi'iniBiaiiaaDBiieiiiiBDDiaaaBi  a»aaaG3iaa  laooaa  a  aaflBBUa  aaSaa  alaa 
nBaaaaa  a  a  ■  a  a  a  a  a  jq  a  a  .:3  a  a  bub  a  a  atiasaSi  a  a  sn^saaa  aUQBa  du  bg  BBa  aSBBU  tisaaa  a  aallaaal  la  •  *  • 
iBlBaalBaaaaBiaaSB'SBBBBBBBaBBBBaaaaaa  BBaBBHsaaaaaa  aBlaGBaao]  fcfioe Blaa aaaBBaaBBBg  ■  •  • 


■  BaiMiaaeaiB  lllii  aalialaalHBiBiillt3i]BBBiii>  ■ 

S^aaBiieaiGo  ^ItlaBHIa  aaa  a  ■■  d  aa  a  ao  nsfliaalaa  a  •  • 
1333  a  a  a  a  bXKIHbIQ?!b&]bJibb  asSsoaaBB  a  eBBnaBllaaBB(  • 


- KflaaiBia  iiZaEQailaalla  aBBliai  tijQUUi£aaBBaj|BBlBaiaBHaaaBlai<Baaiiiaaaoi:BaaBaBaaBBaB.  •  - 

■  aaaaasiilBBlli  a  BilEaaaEaaaaa  a  aH^aaaa  cpsOflaasiBlaa  a  ■  a  BiDSaaBBaiBBiaBDasaUlBalaiaian  a.  • 
•■aaaaai3aalBaBiaaEaEilaiBaaaE|a^aBRaa  aclJBa  USB  I  a  bccibcb  a  as  ao  8  bbQQibBBji  •  a  oa  B2oQ  B  saaCSlai  ■  • 

•  IB . .  a£8la  BBa£}aRUaaaBaRKBBDBDr;S|9BdBlsEilBUaa  tBslana  a  a BansZEiOBalaa  • 

•  iBOBaflBlfl!ri.]BBBEBBB:J£3Bs£?QSj  sBBB  oaQ!iaBB3IlaaBBBBM3Bbi  3Gbna  Baaa  a  BBEtt  8  a  uBBB  a ■  a  .  .  ■  o  n  airn  a  a  B  ■  a  a  a  a  a  • 

■  a  BBaaaaaBBiaa  a  bbbb  b^uaaaaDTlBBkaa  KuainRaBaaESiUBjaDB  JESKiUBBaaaa  aaBBa  maim'  iBiiaiiiaiaaBlii  ■ 

■  •  a  aBBalllaa  cgiaa  8?*7CHagBBaBBaoaa9  aa  a  dUHIa  aaa  cBBSllaa  a  ala . aaaa  a  •  •  •  aaaa.  •  ■  -  •  BBBflBBlflaalBBa  • 

•  •  iiBlBspaaGlflla|j^Bi!t8laaElaal!BlaallBlllaBBiiBBBiaiiaBaBa . . . .  ■  a  a  a  o  sea  a  ti;^aaia  ■  • 

•  «aa  aaa  eUanatiC  as  ar.Bie  Barn  la  BBa  bibbIiIBIbIIbsIIbbI  aaaa . Baa-  .alaa.-  •  .  a  Baa  .  a  a  aBa  aaBailSBlIaa  •  ■ 

•  •  a  B  niiaBIDaa  sr  BaalaBBllFlJBSBBaa  aa  aiacDBaa  uaaaa  a  .  •  .  aCaa •  •  ■  B3.  •  •  ■  aS&Bla  •  •  naBaaiBOBiaaiDUaiUllliii  ■ 

•  •  •  a BBBllaBnaaaQUa  allliO^iBBaa uaalBaaia ■ a ■ . a  a •  • bob* •  -aaa . a B BOB a  a acSaaQBB BaBaoBaaaB alallaa * 

•  a  a  liBBaaaBBeSBBaa  aBBaaQBBaaa  aaBaaaaaa  •  .aaaa.  .  .  .alaaa  -  •  .a  a  Baa  a  ■  aaRBal  BalilBiB  asaBBgflJ  fa  aaiialBli  a  •  ■ 

■  ■BBBaEaaaiaauniaaaaBaaaiaaalaBa.  a  Blaa*  •  •  •  •  aB»  •  •  a  Ilafla  a  iiatiZai  a  a  ana  Baa  all  IB  as  a  aGBaBBRliaiaiaBBaai  ■ 

•  BaBBBBBBaaoaBBSBiaaViaa.  •  •alaa.  •  •  •  a  Dl  a  •  •  • ■ Baa  •  •  •  •  KSBsaBawBBdaEDSBaBaaansReaaQCrlDRBa  lOadEOIaBBRiaai  a  ■ 

. •IBaHiaaiiaaBaBiBaaSBai ■  •  ■ a  aBa  •  -  • aBa  a* . .aBaBiaaiBai 7 aaa aBeS laeBsa BBS  BaaBuBBSssaaaa essseaBliaaBaaai • 

•  •aaiailBi.aBiaia...aiai-  aaaa.  .  . .  aa  a  a  ■  ■  .-9oeBaE§B.  Br.rjriJaogaanET)?  iiaBBNBsaaa  aaa  aBa  a  n  .'.!rr.’JHaaBnan9aaaa  .  ■ 

•  i  a  a  a  a  aaa  •  •  .aDBa  •  •  ■  •  a  a  a  a  «  •  •  .alaa  ■  aa  a  a  •  a  aaa  aloalBEi  aaa  aBBa  a  Ba  alia  a  iSiilaiBBBBaa>aBBaaai  ailiiaiiBiliBiiaaBi  ■ 
•■■■aaaaa*  ••aBBa.  .  .  ■  •aBaiBi.aiaBaiiBBaBBiiHgaaiBaiBBBBailiaaiaaalkBBB|li aaBaaBlDBaaB|liBlBBalBliaiiai a • 

•  a aa aBa a  a •  >aafliM  •  • a  a  a  a o Ba aBa a Baaa BalBa aaB a a CBBa saBlaia aBBa a  Bill (Si. a BBa a aaa a  a a aBBa a aBBBBaall aaaa ■  • 

■  aaBBB  Ba*  •  -•aBaata.ifBaBBaaaBaBEElaaBBBBDaarfllB  .  a  BOB  la  a  la  a  aBBa  •  as;  eh*  aBBa  |  aaa  BBa  a  aBa  BBa*  a  ■  aBalaa  • 

a  a  alQllaa  .  •  ■  BBB  a  aBaBBSiaceSaBaa  calaaBaseBaaELrlBjia  aaSanuaoE53Ba«iaS3B«iB  a  BfiflQlaHBS3qE7!aBB  ■  •  a  ■  .  a  aa aaoonBOBa * 

■  aaos3BaiiiiaaBaEeaB:iBD:aDiiaua]a3aaiBBBaa3aEcJlir.3aalaBiiiiEUDB9iBa3alaiBBBBaiaBia’ljiaB*a  .....  iiatiaaiBlBai 

.  laBinBaaiiDcasaaiiaanRDxaaiiseiicn'jHjiBaaiBs^gajCuiJcRNiBDQSK'saiaaiiiBaiBMii  ■  •  caaa jb aaBBaa 

•  aBaBaBSBBRasusnaaB  BBCi2aauaBB-eaaaaccaBiaaaiBEllaBBaC»Baoa5!lBaaaa  jBaaa  ■ala  . . . . . uaaaa aaaaaa 

- - — - — — - — - - —  —  —  -sSkjlia  J-3 a  a  a  ag Ba  a  c  alia b a  ■  •  a  a  »  o  a  a  ■  ■  ■  •  ■  •  •  .Bllai.  ■  •  •  illBaa  aalBaa 

laBaaaaa. • ■ ■ .8888868881088 
•  ■L’tSiaaBaa  a  ■  •iBaaflliBlaalBB 
• •iHsiaaBaaaaaaBBHSflaalBaaBB 
aBaaaia..  •  .  •  aBecdii  •  ■  ••aaaaia.  •  •  ■  .BltiaaaBaBleDiBaLiBiMalBaBB 

■  =  ?aB.«  •  •  ■  IBM  r  a . .  •  ■■■  ainiaa  laaoiliaaBlaa.  alflaai 

aSlBBBa  .  .  •  ■  a slls ■  . acan^uia  afilBftaDfinaaana  aBBaa  aBsaailBBaa 

■  uaaaa*  •  •■cai a  •  ■  aa  a aBa ansa  Elan  ids a uJ ana a !□■«■ oCH3 a iOIeIIIbIbb 

•  BBa  aa  •  »  •  •  .  .  ■Jilaaa  ■  ibibbb3I)M  ala  Kill.  BTiQD^aiBaaa  sQaaisyaiOUBiBai 

•  aa  a  Baaa  a  ■  aaC^BuBgaa  aBa  a  a  a  arJBaaoBBHSBa^y  aBa  otalliOE  jaffBarjKacaaBaBB 

■  a  .B81lBfflBaaaOiacilB!S93fl.|Diaai3.1CUlGOlBU3BDUD8BaBQOBBtU»aBllElBBI 

Cfls  E j3Sb  a • a  ■  bj|3i  a  ■  "jlBBa"  ■  Ba  aBa  uaHGJaBaaliaatao  a  a  a  aais  BBIBJC3B  BsIBI 
■SaiDlaBB  •«EasgiaBacQHa|BBaHt.BB  aBaaileag oBaaMBaa  Ba  esaa  alsa^JitlBBfliB 


a  Baaa  alBaBa  aBBBaa  aeaGlaalaaa  BsSaoBDonBaai  laalaBaiaiiaiuaflia 

•  a  bi  bbB<]IIIb'l:3BKid  D39Bai'3!]s)9BKdQS3saBla  alia  . . . 

•  BBaaBDiieaiGgBitiMaaBaBBoiaaL'SaaaiBiaaBiBaaaBBBi  •  •  ■  "BBa  .j« 

-  •  ■  BBlBiRcasaa  aU/sa  ji  aBBBaa  Jam  •  •  aaa  aaa  a  ■  ■  •  •  •  air 

■  •aaiis^'.'!jaa.g3vinoiBB(]aiaue;JBDg3M.' 

•  .aBSUOasinsaaBaenDlaiaaMUcSBEla.* 

•  aBBsKBuKaaBaBabdlBBa . . 

■•iRISlBBBaaaiBBBUlBla . a  Biisaa • 

•  •alaaalai..'  •  a  a  afi  Mlifla  •  •  •  •ansas*  • 

.BIIIBIII' '  lacaGBDa.  .aBBatia 

•  a  Biaaa  •  •  ••aiaBai-  •  ■ a alflS ■ a  a > a  a c 

•  ■  CB a  BB a  •  iaSlaaai  .  •  •  aaaana ■ a  a "Boa 


•CBB| 


52! 


_  aSSaaaa..  ..aiiiiBBaaaaBBiBaa^Saia ■ a aaBSifla bSSEbEbsI oBaa  a DUB a a ■SiaaHaEaBaiflBii iia JBBi BBBBaa 

■  BSaiia  ...  laaaaaiaaa  •  i:BaaiE'lfl3aa  aBa«»ao  ■  :^BaBaBB^aaBEetiUB9aanlBaBaaaaacBBaBaaaaia 

.liiaa . alii!  iiiiBaagEiiia  a9|j  cbbb  nfsaE  a  at  a  iSJlliaa  a  f  ^  n  a  n/.la  . .  aa  a  iiiani 

•  a  aaoBaa  a  •  »  a  a  asas  a  B0aB3i  a  a  ta  Ba  ■  a  tin  a  aSlu3  ana  a  b3biu>]  a  a  a  Is  as^aa  ■  a  a  *  ti  EsSsa  a  osa  a  a  a  a  a  a  ■  •  ■■  aBBBa  •  •  •  bed  aBBa  a 

.  i(BBlBaaiiaiajnuicaaiS!nBBaiaBaBBniBCB»aigBa>alaaala>aBiioaaaflDifjaiBaaaiBaisBill«>.  •  •  .  ausia - 

•  a  cIlUBBaa  bDbbUiiESESQIEIdi  a  ■  a  BBiaofljiirDnEsiBBisMUilBBai  aaSSa  sK'|ilB  «  . .  1  ■  ■■■■ 

BaBOIaDBaB|aaaoaBiBfliCZBiaiBQjBS:383l3!iSaaaCl3|UBillEUB.nBaaal3aa  aaa  •  •  ■  mssi.-  •  • 

- "-BBBBBapaatSaolJaaBrS*  BiaisgBBBBBJia  BBiBHaBBBBBB  a  . . .  . . 

_  —  -  nuiijaBaBaBaBBB'*^|aBa  aaaaBOBBBa  •  ■  •  aaajBa.  •  . . 


.  (8888880888 

■aBBaoaQa"' 

- B'.lfll 


■  ■  a  B  a  o  a  '1 


jEEBBarBnaEJoaair, 


i'4anaBaBBBBCB3naaBaag 


mru< 

e£a 

n23 


BOB.  I 


naasc _ 

•  anaai  ■ 

•  •  a  aBa  ■ 


-  -  ■  ailSaBla  a 
.  nBiaaiaa. 
■  a  aSaaaBBa • 
aaSlarnaBa 

IBa  a  •  •  a  aaaBa . aBUtjBBBaaaa-ulB^Bai 

Baa-  .aaaBa.  a  a  a  a  a  au’lcaaaaE  nseiBiaBa  ■ 
aa .  . . . . . 

■  ■  ■  |  i  i  laaBn  sllBSlBaaa  auBKa^BaGga  aa  •  •  •  •  •  .alaa*  •  ■  ■  .aaBSa*  •  •  •  «  •  .llaaBBi  .tBBiaB3BallBaalt:!|'la3laaillB** 

niaaBaaBicSlEaiBBisaBaiai.aBBSBBaa  •  •  ■  •  a  a  ■  a  ■  •  •  •  a  ■  ■  a  B  a  a  .  .  a  ■  a  cIL-SB  a  a  ia;-,ts»i  ifliiscaa  aa  oeTi?  aa  aa  aBBBaa  ■ 

.  . . a  a  a  a  a  .  •  ■  •  acfia.  ■  •  ■iBlliiaiBBBaiBlaiillB3BlianaliJBaaiii  aBBaala  aDBla  a  • 

a  atagilflalBaa  ■  ■  <  .  ■  e  t3  a  a  .  ■  •  •  •  a  a  a  a  a  a  •  ■  ■  •  ■  ca  B  a  ■  ■  iaBXB30iBBaflBBa8IBili|aBBBBaB  S3  BailiilaiiaaaallilSDBa  • 

illilillBlaBa  >  •  •  .a  aBBa  ■  •  •  •  •  a  a  aBBa  •  ■  •  ■  •iiDiBBaalBMBaaallBBBaiacSaaa  IB  aBKaaa  IBBBfla  BB  a  a  ail  a  a  aBBa  ■  • 

.  a  a  laaaaaSl  ■  •  •  . . a  aaa  . . .  alaa  aagflaa  a  a^Baafla  a|aaa  iBai  a  aaaa  BQBBa  a  allaa  aBa  .  • 

.  •  a  a  a  allBlla  •  ■  .a  aBa  •«■••■  bbbbb  a  •  ■  a  a  .i|BaiiBiBalai*lBliaBiBllaBaallBaaBBiaalaal>MaiBBaaaalaa. 

•  ■■! aBlaSSa . •  • . a aaaa a  a ■ a  a aBBBBa a alSBBlfla aa aBBaa aBBBBBa Biaaa a aiaaaa aaiBBaala a aBBl aaaaaaaa aBa a • 

•  ■ aaaaaaaaaa a • • >  aaaaaaaa«BBBaaaB9baJ BllI]aBBBIall|la illaaia|BBlliaiilllBBaiiiaalllillaaiailai . 

■  laBaaaaaBBa* •aalaaaaalaaBaaaa  aBlaa  a laalaiRllaiBBBa alBBBliaa l ■ bbb a  a  aBa ■ • * ' 

.  .  iiiiiiitDlslaaaBBBaBlBlsaBaBBaii  99  eiaBBaaBiElaaBBaai  aaaaaaaa  aa  IBBaaEBaaalBBBBaBiaiqEaviiSiaii 
•  ■■■iBiaasiiaaiiBiSBiBBBBkCEisiyusQSansaEnQE^’iBflaiiaiBBiaa  ■aDiBBailBBaaBaiBla  aaaalaaiRailBa  •  • 

.  •BBiSBi?.iaBiaBiiB(95laiBBai!S,9alialcsiatiai3lliacilaBl|aBaaiBiill8fiiaaaaaiiaaBiiBaii;:seilai' 
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inaBBHBDIaaaBEiBaa 
aaaaa aa BIB Baaaa aaa 
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!■■■■■■■! 


HOLE  DIAMETERS  (INCHES)  ARE  (TOP  TO  BOTTOM)  1/32,  3/32,  5/32, 

1/8  AND  1/16. 


FIGURE  32.  HIGH  SPEED  RESOLUTION 
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FIGURE  33.  MEDIUM  SPEED  RESOLUTION 
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FIGURE  34.  LOW  SPEED  RESOLUTION 
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B.  Precision 

The  critical  factor  for  the  overall  precision  of  the 
system  is  not  so  much  the  precision  associated  with  the 
machine  as  with  the  ability  to  remeasure  the  patient  at 
exactly  the  same  site  since  both  the  amount  of  trabecular 
bone  and  its  density  vary  significantly  at  sites  axially 
along  the  arms  or  legs.  This  means  that  precision  of  the 
system  is  predominantly  a  function  of  the  ability  to 
determine  the  site  of  measurement.  To  reduce  these  errors  in 
precision  due  to  repositioning  use  of  the  scout  scan  was 
i nst i tuted . 

A  number  of  repeat  scans  were  made  on  a  40  year  old 
male  normal  (52).  Each  scan  involved  repositioning  the  arm 
using  the  scout  scan  and  making  three  measurements  on  the 
arm.  The  first  measurement  was  made  at  the  site  selected  by 
the  scout  scan  while  the  other  two  measurements  were  made  at 
sites  two  millimeters  of  either  side.  These  three 
measurements  were  then  interpolated  to  produce  the  same  area 
of  bone  for  analysis  purposes  as  was  used  in  prior  analyses. 
The  data  obtained  from  these  measurements  (Table  2)  shows 
that  a  precision  of  better  than  21  %  is  achievable  for  both 
cortical  and  trabecular  bone  density. 

C.  Accuracy 

Accuracy  of  the  scanner  is  a  difficult  quantity  to 
define  since  it  relates  to  a  "correct”  value  being  obtained. 
For  measurements  of  bone  this  means  that  the  bone  would  have 


. 
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PRECISION 

EVALUATION 

Measurement 

Scan  * 

Bone  Area 

Trabecular  Bone 

Cortical 

Number 

Number 

(Pixels) 

Density 

Densi ty 

1 

1 

3362 

0.503 

1.301 

2 

1 

3577 

0.491 

1.250 

2 

3293 

0.519 

1.333 

3 

1 

3552 

0.501 

1.250 

2 

3248 

0.522 

1.328 

4 

1 

3389 

0.509 

1.275 

2 

3155 

0.  535 

1.370 

5 

1 

3586 

0.497 

1.237 

2 

3308 

0.510 

1.330 

With  the 

above  scans 

i nterpol ated 

so  as  to  produce  a 

bone  area 

of  3362  the  results  are  as  follows 


1 

3362 

0.503 

1 .301 

2 

3362 

0.512 

1.313 

3 

3362 

0.514 

1.299 

4 

3362 

0.512 

1.286 

5 

3362 

0.507 

1.292 

These  interpolated 

results,  when 

considering 

a  normal  distribution. 

result  in  standard 

deviations  of 

0.884  7o  for 

trabecular  bone  density 

and  0.847  %  for  cortical  bone  density. 


*Note:  Although  three  scans  were  made  during  each  measurement 
only  those  two  scans  that  were  used  in  the  interpolation  are  included 

TABLE  2. 
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to  be  removed  and  "ashed"  to  determine  the  true  mineral 
content.  Obviously,  this  is  only  possible  to  perform  with 
cadavers  and  consequently  the  sources  of  data  are  limited. 

To  circumvent  these  problems  a  solution  of  di -Potassium 
Hydrogen  Phosphate  (K2HPO4)  was  used  (53).  The  solution  was 
mixed  to  a  Known  concentration  and  placed  inside  a 
plexiglass  cylinder  (Figure  35).  CT  measurements  were  made 
on  a  number  of  different  concentrations  and  the  resultant 
linear  attenuation  coefficients  of  the  image  were  compared 
to  the  Known  values  of  the  dilute  solution.  Results  of  these 
tests  are  shown  in  Table  3  and  Figure  36.  These  results  show 
that  although  the  absolute  values  may  not  be  entirely 
accurate,  the  percentage  change  in  value  is  quite  accurate. 
Since  the  mineral  content  present  in  bone  does  vary 
tremendously  from  person  to  person  it  is  measures  of 
percentage  changes  that  are  of  clinical  importance  rather 
than  absolute  values  of  mineral  density.  One  reason  that  the 
absolute  values  are  of  limited  accuracy  is  that  the  linear 
attenuation  coefficients  for  the  solution  are  calculated  on 
the  basis  of  a  single  energy  which  does  not  coincide  exactly 
with  the  energy  spectrum  of  Iodine  125. 

D.  Radiation  Dose 

Calculations  of  a  theoretical  radiological  dose  (as 
performed  on  the  following  pages)  show  that  the  total 
radiation  dose  received  during  the  course  of  a  scout  scan 
and  three  adjacent  scans  is  0.9  rads. 


' 
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MEASUREMENT 
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FIGURE  35.  TEST  CYLINDER 
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ACCURACY  EVALUATION 


Measured  Linear 

Attenuation 

Coefficient 

( per  cm. ) 

Percent  of 

Maximum  Measured 

Linear  Attenuation 

Coefficient 

Calculated 

Linear  Atten. 

Coefficient 

( per  cm. ) 

Percent  of 

Maximum  Calculated 

Linear  Attenaution 

Coeff ici ent 

1.882 

100 

1.709 

100 

1.852 

98.4 

1.675 

98.0 

1.814 

96.4 

1.641 

96.0 

1.776 

94.4 

1.606 

94.0 

1.765 

93.8 

1.598 

93.5 

1.747 

92.8 

1.589 

93.0 

1.739 

92.4 

1.581 

92.5 

1.735 

92.2 

1.572 

92.0 

1.690 

89.8 

1.538 

90.0 

1.654 

87.9 

1.504 

88.0 

1.615 

85.8 

1.470 

86.0 

A  linear  regression  of  the  percent  calculated  linear  attenuation 
coefficient  was  run  against  the  percent  measured  linear  attenuation 
cefficient  (using  the  SPSS  package)  and  found  to  have  a 
correlation  of  r=0.9987  and  a  significance  of  p  <  0.0001 


TABLE  3. 
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FIGURE  36.  SYSTEM  ACCURACY 
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Radiological  Dose 

Assume  that  an  arm  of  diameter  60  mm  is  placed  in  the 
center  of  the  scan  area.  This  results  in  a  distance 
of  43  mm  from  the  source  to  the  nearest  point  on  the  arm. 


At  a  distance  of  43  mm  the  source-collimator  configuration 
was  found  to  produce  the  following  approximate  intensity  distribution. 


The  scanner  moves  a  distance  of  106.2  mm  in  12.89  sec 
resulting  in  an  average  speed  of  8.24  mm  /  sec.  The  area 
irradiated  at  100  %  intensity  consists  of  a  length  L,  where 
L  =  2  tan  ( 23 °/2)  (43  mm)  =  17.5  mm.  The  area  irradiated 
at  an  average  intensity  of  1/2(100%  +  50%)  =  75%  consists 
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of  2  lengths,  each  one  of  length  LI,  where 

LI  =  1/2  (2  tan  ( 60 °/2)  43  nm  -  17.5  mm)  =  16.08  mm. 

Thus  the  total  length  irradiated  at  effectively  100°/  is 
17.5  (100%)  +  2(  1 6. 08) (75%)  =  41.61  mm  which  at  a  speed  of 
8.24  mm  /  sec  means  an  effective  (100%)  irradiation  time  of 
41.61  mm  /  8.24  mm/sec  =  5.05  seconds. 

This  is  the  irradiation  time  for  one  translation.  However, 
each  scan  consists  of  7  translations  made  at  varying  angles. 


The  rotation  between  angles  is  25.71°  and  the  approximate  hotspot 
(area  subjected  to  maximum  radiation)  is  that  point  on  the  arm 
perpendicular  to  translation  number  4  (denoted  by  point  x) . 

Thus  at  y=  90°the  irradiation  is  100%  while  at  y=  0°the  irradiation 
is  0  %.  Hence  irradiation  is  proportional  to  sin  y. 


Transl ation 

1 

y 

=  90°  -  3(25.71°) 

=  12.87°  sin 

y  =  0.22 

Translation 

2 

Y 

=  90°  -  2(25.71°) 

=  38.58°  sin 

y  =  0.62 

Translation 

3 

y 

=  90"  -  1(25.71°) 

=  64.29°  sin 

y  =  0.90 

Translation 

4 

y 

=  90° 

sin 

o 

o 

n 

>- 

Translation 

5 

y 

=  64.29° 

si  n 

y  =  0.90 

Translation 

6 

y 

=  38.58° 

si  n 

y  =  0.62 

Translation 

7 

y 

=  12.87° 

sin 

y  =  0.22 

4.48 


Therefore,  for  one  scan  the  total  irradiation  time  (100%  intensity) 
at  the  hotspot  is  4.48  (5.05  sec)  =  22.62  seconds. 
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A  new  source  of  1.5  Ci  I125  was  measured  using  an 
ionization  chamber  at  a  distance  of  0.5  meters  (so  that  the 
entire  chamber  was  irradiated)  and  found  to  have  a  dose  rate  of 
11  mr/min.  Using  the  inverse  square  law  the  dose  rate  at  43  mm 
would  be  11  mr/min  (500/43) 2  =  1487  mr/min  =  24.79  mr/sec. 

This  means  that  the  total  dose  per  scan  at  the  hotspot  is 

24.79  mr/sec  (22.62  sec)  =  560  mr.  Since  the  speed  of  the  scout 

scan  is  the  same  as  for  a  normal  scan  the  effective  (100%) 

irradiation  time  is  5.05  sec  leading  to  a  dose  of 

5.05  sec  (24.79  mr/sec)  =  125  mr.  Assuming  adjacent  measurements 

produce  a  50  %  increase  in  dose  the  total  exposure  dose  for 

1  scout  scan  and  3  adjacent  scans  is 

1.5(560  +  125)  =  1027.5  mr  *  1  r  (roentgen). 

At  30  KeV  this  is  approximately  equal  to  an  absorbed  dose  of  0.9  rad. 
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Although  these  theoretical  calculations  do  provide  a 
rough  estimate  of  dose  there  are  a  number  of  associated 
problems.  One  major  problem  is  the  assumption  that 
absorption  of  radiation  is  only  significant  at  the  surface. 
Although  this  may  be  true  for  denser  objects  with  high 
linear  attenuation  coefficients  it  is  not  true  for  skin.  The 
linear  attenuation  coefficent  of  skin  is  approximately 
0.45/cm  (for  30  KeV )  which  means  that  some  1 . 5  cm  of  skin 
must  be  traversed  before  the  intensity  of  the  radiation  is 
down  by  a  factor  of  2.  Since  there  is  a  limited  amount  of 
skin  at  the  forearm,  radiation  passing  through  skin  only 
will  have  a  significant  effect  on  the  dose  rate  whether  or 
not  it  directly  impinges  upon  the  area  of  consideration. 
Although  most  rays  do  pass  through  some  bone  and  thereafter 
do  not  contribute  appreciably  to  skin  dose  it  should  be 
realized  that  there  is  an  element  of  variation  from  person 

to  person. 

Another  area  of  concern  is  the  effects  of  scattered 
radiation  from  the  bones  and  their  effect  on  the  dose  at  a 
particular  point.  It  is  relatively  complicated  to  account 
for  scatter  on  a  theoretical  basis  and  as  such  the  use  of 

t  Detectors  (TLD's)  has  been  suggested. 


Thermo- L urn i nescen 
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Thermo- Lumi nescent  Detectors 

TLD' s  are  capsules  of  powder  or  crystals  that  react 
chemically  when  exposed  to  radiation.  After  exposure  to 
radiation  the  powder  is  removed  from  the  capsules  and  heated 
up.  When  heated  the  powder  produces  photon  flashes  the 
number  and  intensity  of  which  are  proportional  to  the 
radiation  exposure. 

It  is  presently  being  proposed  that  a  number  of  holes 
be  drilled  in  a  plexiglass  phantom  and  filled  with  TLD 
powder.  Subsequent  to  making  a  number  of  measurements  on  the 
phantom  the  powder  is  to  be  removed  and  tested  to  determine 
the  radiation  exposure.  The  plexiglass  phantom  will  contain 
aluminum  tubes  to  simulate  the  radius  and  ulna  and  will  be 
made  in  the  approximate  size  and  shape  of  a  human  forearm. 
Measurements  of  dose  performed  in  this  manner  should  produce 
a  profile  of  the  radiation  dose  throughout  a  human  forearm. 

E.  Contrast  Resolution 

Contrast  resolution  could  be  measured  by  examining  the 
limiting  spatial  resolution  with  differing  contrast  media  in 
the  holes.  One  possible  contrast  media  is  di-Potassium 
Hydrogen  Phosphate  (K^HPOy)  mixed  in  a  water  base.  If 
solutions  of  varying  concentration  were  placed  in  the 
1/16  th  of  an  inch  diameter  holes  then  the  system  limited 
contrast  resolution  would  be  defined  by  comparing  the 
smallest  linear  attenuation  coefficient  of  the  imageable 
solutions  with  that  of  plexiglass.  The  major  difficulties  in 
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these  measurements  will  be  placing  the  proper  concentration 
solutions  in  the  holes  and  ensuring  that  no  air  bubbles  are 
present.  Measurements  of  contrast  resolution  have  not  been 
made  due  to  these  difficulties  and  the  limited  usefulness  of 
such  measurements. 


VI.  System  Limitations  and  Future  Plans 
Four  of  the  major  limitations  at  present  are:  (1) 
statistical  limitations  due  to  the  use  of  a  radioisotope 
photon  source,  (2)  a  lengthy  scanning  time  due  to  use  of  a 
second  (versus  a  third  or  fourth)  generation  CT  scanner ,  (3) 
a  lack  of  flexibility  in  regards  to  scanner  control ,  and  (4) 
an  inability  to  view  the  image  immediately  after  scanning 
due  to  limited  processing  facilities.  This  chapter  will 
explore  these  limitations  and  will  suggest  ways  to  reduce  or 
eliminate  them. 

A.  X-Ray  Tube 

Subject  to  appropriate  funding  present  plans  include 
the  purchase  of  a  mini-X-ray  tube.  As  a  result  of  physical 
considerations  photon  flux  is  more  limited  in  the  case  of  a 
radioisotope  source  than  for  an  X-ray  tube.  Although  use  of 
an  X-ray  tube  may  result  in  higher  radiation  dose  rates  it 
will  also  facilitate  a  reduced  scan  time  and  hence  the  net 
radiation  exposure  may  not  be  radically  altered.  With  the 
aid  of  filters  the  X-ray  tube  will  also  allow  a  greater 
selection  of  photon  energies  to  be  used.  This  may  allow  the 
scanner  to  perform  elemental  identifications  through  the  use 
of  differential  energy  scans  and  to  be  used  for  scans  of 
areas  of  high  radiological  attenuation  such  as  the  knee.  The 
use  of  filters  may  also  reduce  the  problems  of  beam 
hardening  caused  by  the  multi-energy  X-ray  radiation. 
Although  beam  hardening  will  be  a  problem,  so  long  as  the 
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"raw"  data  is  available,  it  is  possible  to  correct  for. 

B.  System  Three 

Present  plans  include  the  construction  of  a  generation 
three  system.  This  improvement  will  result  in  a  substantial 
decrease  in  scan  time  hence  reducing  movement  artifacts, 
improving  patient  comfort  and  reducing  the  dose  of 
radiation. 

C.  Control  System 

The  major  constraint  at  present  is  that  the  control 
programs  are  " locked "  into  EPROMs.  In  order  to  alter  any  of 
the  scan  parameters,  such  as  speed,  scan  length,  distance 
between  projections  and/or  the  angle  between  profi ies,  a  new 
program  must  be  written  into  the  EPROMs,  o  circumvent  this 
problem  it  is  proposed  that  the  EPROMs  be  programmed  with  a 
number  of  basic  subroutines  and  that  all  control  parameters 
be  loaded  into  the  MPU  from  external  sources.  One  of  the 
subroutines  contained  in  the  EPROMs  would  be  a  program  that 
reads  data  into  the  MPU  via  the  PIA  (M6820).  Data  to  be  read 
into  the  MPU  would  be  loaded  into  the  HP  2100  from  either 
the  console  or  a  series  of  programs  contained  on  magnetic 
tape.  This  data  will  be  subsequently  transferred  from  the  HP 
2100  into  a  memory  module  contained  in  the  CAMAC  system.  The 
transfer  of  data  from  the  HP  to  the  memory  module  (4000  16 
bit  words)  will  free  the  HP  to  perform  other  functions 
subsequent  to  loading  the  CAMAC.  Control  parameters  will  be 
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loaded  from  the  memory  module  to  the  512  bytes  of  RAM 
(M6810)  contained  in  the  MPU.  Direct  interrupts  will  also  be 
provided  between  the  MPU  and  the  HP  2100  computer .  The  read 
program  which  is  the  subroutine  that  reads  data  from  the 
memory  module  to  the  mi croprocessor  unit  has  been  assembled 
and  is  included  in  the  appendix. 

D.  Processing  Facilities 

A  new  32  bit  minicomputer  has  been  ordered  (SEL  32/57) 
and  delivery  is  expected  for  the  fall  of  1980.  Purchase  of 
an  array  processor  dedicated  to  the  reconstruct  ion  is  also 
being  considered  and  appears  likely  at  this  point  in  time. 
The  resultant  upgrading  of  computing  facilites  will  allow 
the  image  to  be  displayed  on  the  Grinnell  imaging  system  at 
the  same  time  or  shortly  after  a  scan  is  performed.  This 
means  that  the  validity  of  the  scan  can  be  checked 
immediately  and  should  another  scan  be  required  the  patient 


will  still  be  present . 
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125  125 

I  decays  via  electron  capture  to  Te  and  in  doing 
53  52 

so  gives  off  a  35.4  KeV  gamma  ray.  This  gamma  ray  will: 

1)  appear  on  the  outside  .0666  times  per  disintegration; 

2)  cause  a  k  alpha-1  x-ray  .7615  times  per  disintegration  (27.4  KeV); 

3)  cause  a  k  alpha-2  x-ray  .3906  times  per  disintegration  (27.2  KeV); 

4)  cause  a  k  beta-1  x-ray  .2056  times  per  disintegration  (30.9  KeV); 

5)  cause  a  k  beta-2  x-ray  .0426  times  per  disintegration  (31.8  KeV); 

6)  cause  an  L  x-ray  .2226  times  per  disintegration  (3.7  KeV). 

Thus  in  general  the  original  gamma  ray  energizes  an 
electron  in  the  k  shell  causing  it  to  move  out  of  its  orbit. 

This  gap  is  filled  by  a  loose  electron  which,  in  falling  into 
the  k  shell,  gives  off  an  x-ray  of  approximately  27.5  KeV. 

The  associated  half-life  is  60.2  days. 
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APPENDIX  B. 


CONTAINS  ALL  MAIN  PROGRAMS  AND  ALL  SUBROUTINES  (EXCEPT  READ) 
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